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PREFACE 


This  report  was  prepared  by  Hughes  Research  Laboratories, 

Malibu,  California  under  Contract  F30602-76-C-0021.  It  summarizes 
the  results  of  a six-month  experimental  and  analytical  study  of  the  inter- 
action between  multidither  COAT  systems  and  target-produced  amplitude 
modulations  in  a COAT  receiver.  An  interim  report  (Dec.  1975)  describes 
the  early  contract  progress  in  detail.  This  report  deals  mainly  with 
results  obtained  during  the  last  half  of  the  contract,  while  summarizing 
the  entire  contract.  The  principal  scientist  is  Dr.  James  E.  Pearson 
and  the  principal  investigators  are  Dr.  Melvin  E.  Pedinoff  and  Mr.  Stan 
A.  Kokorowski.  This  project  is  a part  of  the  Adaptive  Optics  program  in 
the  Opto-Electronics  Department,  managed  by  Dr.  Viktor  Evtuhov,  at  the 
Hughes  Research  Laboratories. 
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SUMMARY 


The  objectives  of  this  program  are  threefold:  (1)  to  determine  the 

characteristics  of  receiver  modulations  produced  by  backscatter  from 
realistic  targets  that  interfere  with  the  proper  operation  of  a multidither 
COAT  system  and  to  assess  the  degradation  they  produce  in  COAT  system 
convergence  and  performance;  (2)  to  develop  a unified  theory  of  the  inter- 
action between  a multidither  COAT  system  and  target-speckle  induced 
receiver  modulations;  (3)  to  determine  an  optimum  COAT  system  design 
that  is  insensitive  to  speckle-induced  receiver  modulations  and  to  define 
the  critical  design  parameters.  This  report  summarizes  the  entire  con- 
tract work  performed  from  1 August  1975  to  1 February  1976,  concentrat- 
ing mainly  on  the  results  obtained  from  ?1  October  1975  to  1 February  1976. 
We  have  experimentally  studied  the  speckle  patterns  generated  by 

scaled,  realistic  target  models  consisting  of  metallized  or  painted  spheres, 

0 

cylinders,  and  cones  using  a helium-neon  (6328  A)  laser  and  a target- 

speckle  measurement  apparatus  specifically  assembled  for  this  program. 

Target  signatures  were  recorded  with  two  receiver  aperture  diameters, 

one  equal  to  the  transmitter  diameter  and  one  equal  to  one-third  of  the 

transmitter  diameter.  The  larger  aperture  size  significantly  reduced  the 

observed  speckle  contrast  ratio  and  reduced  the  amplitude  of  the  higher 

speckle  spatial  frequency  components. 

Speckle  pattern  data  have  been  reduced  using  a fast-Fourier- 

transform  to  obtain  spatial  frequency  spectra.  The  same  spectrum  is 

obtained  by  translating  the  receiver  as  by  rotating  the  target.  Target 

speckle  signatures  taken  with  most  of  the  metallized  target  shapes  are 
. . . >!< 
very  similar.  The  maximum  contrast  ratio  observed  in  all  of  the 

experiments  was  0.  89  with  the  most  commonly  observed  maximum 

contrast  ratio  being  0.  56. 


Defined  as  (I  - I . )/ (I  + I ),  where  I = observed  irradiance. 
max  min  max  mm 
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Real-time  spectrum  analyzer  spectra  obtained  by  rapid  rotation 
of  the  targets  are  identical  to  spatial  frequency  spectra  obtained  by  per- 
forming a fast-Fourier-transform  (FFT)  on  the  speckle  pattern,  data. 
Increasing  the  target  rotation  speed  caused  the  speckle  spectrum  to 
broaden  toward  high  frequencies  in  agreement  with  the  FFT  spatial  fre- 
quency results.  This  result  is  useful  for  explaining  the  effect  of  target 
rotation  on  COAT  system,  convergence  level.  The  speckle  signatures 
produced  by  rotation  of  the  metallized  targets  are  very  similar  to  those 
produced  on  an  earlier  program  {Contract  F30602-75-C-0001)  by  rotating 
scotchlite  targets. 

Very  little  difference  in  speckle  properties  was  observed  with 
spherical,  conical,  or  cylindrical  targets.  Consequently,  the  metallized 
spherical  target  was  used  as  a vehicle  to  test  COAT  system  performance 
because  it  was  the  most  convenient  experimentally. 

Following  the  detailed  measurements  on  speckle  patterns,  the 
spherical  target  was  used  as  a target  in  the  DARPA/RADC  18-channel 
COAT  system.  The  system  convergence  behavior  was  measured  as  a 
function  of  target  rotation  speed  for  several  target  orientations.  Signifi- 
cant degradation  of  the  converged  power  to  a normalized  level  of  P = 0.  31 
was  observed  when  the  target  rotated  at  = 2ir  rad/ sec  about  an  axis  at 
90°  to  the  incident  beam  (P  = 1.0  at  ft  = 0).  Slower  or  faster  rotation 
rates  produced  less  degradation.  When  the  COAT  system  feedback  loop 
gain  was  reduced  by  15  dB  (rotation  optimization),  the  converged  power 
level  was  0.  93.  Similar  experiments  carried  out  with  the  target  rotation 
axis  at  45°  to  the  incident  beam  yielded  strong  degradation  (P  = 0.  50)  and 
good  recovery  (P  = 0.  80)  of  COAT  convergence  for  the  high  loop  gain 
(static  optimization)  and  low  loop  gain  (rotation  optimization)  cases, 
respectively.  When  the  target  axis  was  parallel  to  the  laser  axis  (0°  ), 
negligible  degradation  of  system  convergence  was  observed,  even  with 
high  loop  gain, 

The  convergence  time  of  the  system  was,  of  course,  affected  by 
the  loop  gain.  In  the  presence  of  speckle  modulation,  the  COAT  system 
oscillates  between  two  levels  whose  mean  decreases  as  the  speckle 
modulation  power  in  the  dither  band  increases.  The  convergence  time 


to  this  mean  level  for  the  condition  of  high  loop  gain  is  essentially 
unaffected  by  the  modulations.  When  the  system  loop  gain  is  reduced 
(AG  = -15  dB)  to  counteract  the  effects  of  speckle  modulation,  however,  the 
convergence  time  increases  at  most  by  a factor  of  6 and  on  the  average  by  a 
factor  of  4.  (The  product  of  loop  gain  and  convergence  time  is  approxi- 
mately constant.  ) 

Although  most  of  the  experimental  measurements  were  made  with 
no  propagation  distortions,  we  have  also  investigated  the  ability  to  per- 
form turbulence  and  blooming  compensation  in  the  presence  of  speckle. 

The  distortions  existed  only  on  the  outgoing  beam  path,  however,  and 
thus  did  not  affect  the  speckle  itself  in  any  way.  In  general,  we  find  that 
the  existence  of  speckle  modulations  does  not  significantly  remove  the 
COAT  system's  ability  to  correct  for  blooming  and  turbulence,  although 
the  amount  of  correction  and  the  peak  target  irradiance  are  usually  both 
less  when  speckle  effects  are  present.  The  degree  of  reduction  in  com- 
pensation depends  on  the  amount  of  speckle-induced  power  within  the  COAT 
system  dither  band  as  might  be  expected.  In  each  case  studied,  the  COAT 
system  could  not  improve  the  beam  strehl  ratio  over  that  set  by  the  speckle. 
For  example,  if  speckle  alone  produces  a strehl  ratio  of  0.  3,  and  turbu- 
lence alone  (when  uncorrected)  produces  a strehl  ratio  of  0.4,  the  strehl 
ratio  of  the  beam  with  both  effects  is  0.  3.-  If  the  speckle  strehl  ratio  is 
0.  6,  however,  then  the  turbulence-plus -speckle  ratio  will  be  somewhere 
between  0.4  and  0.  6,  usually  closer  to  the  0.  6 level  as  set  by  the  speckle. 

Analytical  work  conducted  throughout  this  program  consisted  in  part 
of  computer  simulation  studies  of  COAT  system  convergence  in  the  presence 
of  broadband  speckle-like  noise.  Computer  simulations  of  multidither  COAT 
performance  have  been  made  using  both  an  ad  hoc  "speckle"  spectrum  and 
theoretical  target-  speckle  data  generated  by  the  General  Research  Corpora- 
tion (GRC).  The  ad  hoc  spectrum  consisted  of  23  discrete  frequencies  of 
random  amplitude  and  phase  placed  within  the  servo  bandwidth  of  each  of  18 
COAT  control  channels.  The  data  supplied  by  GRC  were  calculated  using  a 
full  diffraction  calculation  and  the  following  parameters:  target-range  = 2 km, 
shape  = 10  cm  diameter  segment  of  1 m radius  sphere,  illumination  = uniform, 
material  = pyroceram,  assumed  completely  random  in  phase  reflection;  other, 
wavelength  = 10.6  pm,  annular  receiver  (1.2  m x 1.0  m diameters)  about  a 
20  cm  diameter  transmitter. 
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Both  types  of  modulation,  applied  as  multiplicative  receiver 
lio'.'je,  gave  qualitatively  similar  results.  That  is,  as  the  noise  power 
increased,  the  average  convergence  level  decreases  and  the  COAT  system 
convergence  level  (strehl  ratio)  oscillates  between  two  well-defined  levels. 
Qualitatively,  this  observation  also  agrees  with  experiment.  Three  cases 
were  run  with  GRC-supplied  data  corresponding  to  target  rotations  about 
an  axis  perpendicular  to  the  illumination  beam  direction  of  £2  = 0.4,  2, 
and  10  rad/sec.  The  average  convergence  levels  from  these  runs  were 
0.  74,  0.  32,  and  0.  51,  compared  to  0.  95  for  no  speckle  modulations. 

An  analysis  has  been  developed  to  aid  in  understanding  the  experi- 
mental and  simulation  data  and  to  provide  a means  for  predicting  COAT 
system  performance  :in  the  presence  of  speckle  modulations.  The  analy- 
sis includes  a statistical  model  of  a multidither  COAT  system  and  incor- 
porates speckle  modulations  as  a noise  °ource.  The  analysis  defines  the 
system  performance  in  terms  of  two  parameters:  Cg,  called  the  "speckle 
coefficient"  and  p ' , the  ratio  of  the  average  dither  signal  power  to  the 
average  speckle  modulation  power  after  synchronous  detection.  The 
quantity  Cg  is  proportional  to  the  square-root  of  the  total  speckle - 
produced  power  within  a servo  bandwidth  around  all  the  dither  frequencies 
and  can  have  a maximum  value  of  \fz' « 0.  707. 

The  analysis  predicts  that  the  observed  maximum  and  minimum 
convergence  levels  are  defined  by  p'  = 0.  5 and  p'  = 2.  0,  with  p'  = 1.  0 
being  the  average.  It  also  says  that  the  £2  = 2.  0 rad/sec  case  supplied 
by  GRC  does  not  represent  the  theoretical  worst-case  since  Cg  = 0.  24 
for  that  scenario  (although  given  the  COAT  parameters,  it  is  the  worst 
case  of  the  ones  studied).  Very  good  qualitative  agreement  has  been 
obtained  between  the  analysis  and  computer  simulation  data  and  with 
experimental  data.  Quantitatively,  the  experimental  and  simulation  data 
agree  when  the  parameters  are  scaled  properly.  The  analysis,  however, 
is  somewhat  pessimistic  in  its  predictions.  For  example,  the  £2  = 

2 rad/sec  GRC  case  most  closely  scales  to  a laboratory  experiment  where 
£2  = 2tt  rad/sec.  The  experimental  strehl  ratio  is  0.31,  the  computer 
simulation  result  is  0.30,  and  the  analysis  predicts  0.  18. 


This  program  has  thus  established  both  experimentally  and  by 
computer  simulation  that  speckle  effects  can  significantly  degrade  the 
performance  oc  a multidither  COAT  system  and  we  have  developed 
analytical  and  computer  simulation  tools  that  enable  us  to  predict  COAT 
system  performance  beyond  the  particular  18-channel  system  studied  here. 
We  have  also  established,  however,  that  the  degradation  produced  by 
speckle  depends  strongly  on  the  speckle  contrast  ratio  and  on  the  rate  at 
which  the  speckle  pattern  moves  by  the  COAT  receiver.  These  parameters 
in  turn  are  sensitive  functions  of  the  target  shape,  size,  surface  rough- 
ness and  motion  as  well  as  of  the  target  range  and  the  illumination  laser 
bandwidth.  The  COAT  system  parameters  also  play  an  important  role. 

For  example,  we  have  noted  that  lower  gain  (lower  servo  bandwidth) 
reduces  speckle  effects;  we  have  also  established  that  a higher  dither 
modulation  index  (up  to  a point)  and  the  presence  of  a hard  limiter  (clipper) 
in  the  servo  improve  performance  in  the  presence  of  speckle. 

Additional  work  still  needs  to  be  done  to  quantitatively  establish 
the  effect  of  modulation  index,  of  the  number  of  channels,  and  of  allowing 
the  speckle  itself  to  change  as  the  COAT  convergence  level  changes.  The 
effect  of  servo  bandwidth  variation  with  convergence  level  has  not  been 
included  in  the  analysis  and  such  an  addition  would  increase  its  usefulness. 
Perhaps  most  important  is  the  need  to  establish  what  scenarios  can  pro- 
duce significant  speckle  modulations  that  fall  with  the  COAT  system  dither 
band  and  how  probable  is  the  occurrence  of  such  scenarios.  If  the  prob- 
ability of  occurrence  is  high,  then  advanced  techniques  for  reducing  the 
degrading  effects  of  speckle  should  be  pursued. 
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INTRODUCTION 


A.  PROGRAM  OBJECTIVES 

The  overall  purpose  of  this  program  is  to  produce  a consistent, 
unified  theory  on  multidither  COAT/ target  interactive  effects  that  correlate 
experimental  data  with  analytical  models.  Specifically  there  are  three 
primary  objectives.  The  first  objective  is  to  determine  what  kind  of 
modulations  interfere  with  the  proper  operation  of  a multidither  COAT 
system.  The  significant  parameters  of  the  modulation  are  contrast  ratio 
and  frequency  spectrum.  The  second  objective  is  to  determine  the  char- 
acteristics of  receiver  modulations  produced  by  backscatter  from  real 
targets  in  expected  operational  scenarios;  that  is,  to  determine  which 
scenarios  can  produce  the  modulations  which  will  affect  the  COAT  system. 
The  third  objective  is  to  determine  an  optimum  COAT  system  design  that 
is  insensitive  to  speckle-induced  receiver  modulations  and  to  define  the 
critical  design  parameters. 

B.  RESEARCH  PROGRAM  PLAN 

The  18-element  multidither  COAT  system  developed  and  tested  on 
DARPA/RADC  Contract  F30602-73-C-0248  is  used  for  laboratory  experi- 
ments on  this  program,  which  is  an  in-depth  investigation  of  speckle 

phenomena  in  multidither  COAT  systems. 

The  research  program  for  this  contract,  illustrated  in  Table  1, 
runs  from  1 August  1975  to  1 February  1976.  Because  of  ongoing  measure- 
ments on  the  COAT  Measurements  and  Analysis  program  (Contract  F30602- 
75-C-0001),  a separate  target-speckle  measurement  apparatus  was  used. 
The  apparatus  included  a helium-neon  6328  A laser  for  target  illumination 
and  a nearly  collinear  receiver  system  for  target  backscatter  character- 
ization studies.  Theoretical:  speckle  simulation  studies  were  also  conducted 


TABLE  1.  COAT/TARGET-SIGNATURE  RESEARCH 
PROGRAM  PLAN 
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Some  of  the  computer  simulation  work  used  target  speckle-signatures 
generated  by  General.  Research  Corporation  (GRC)  under  Contract 
F3602-75-C-0086.  The  program  is  divided  into  three  tasks:  Target 

Backscatter  Measurements,  COAT  System  Measurements,  and  Analytical 
Studies.  Task  1 represents  the  acquisition  and  characterization  of  targets 
and  target  speckle  data  not  previously  available  to  the  program.  Task  2 
determines  COAT  system  performance  in  the  presence  of  turbulence, 
blooming,  and  speckle  noise  modulation  from  these  targets.  Task  3 is  a 
computer  simulation  study  designed  to  corroborate  or  predict  the  experi- 
mental findings.  Preliminary  results  from  Tasks  1 and  3 are  reported  in 
the  Interim  Report  on  this  contract,  dated  December  1975. 


C.  REPORT  ORGANIZATION 


This  report  serves  two  functions:  (1)  it  presents  the  accomplishments 

of  the  second  half  of  this  program,  and  (2)  it  summarizes  the  principal 
accomplishments  of  the  entire  contract.  The  report  is  thus  both  a second 

half  technical  report  and  a final  technical  report. 

Section  II  reviews  the  target  speckle  characterization  accomplish- 
ments obtained  during  the  first  half  of  this  program  and  presents  complete 
data  on  the  performance  of  COAT  systems  in  the  presence  of  speckle;  bloom- 
ing and  speckle;  and  blooming,  turbulence,  and  speckle. 

Section  III  summarizes  the  results  of  the  analytical  efforts  directed 
toward  investigating  complex  target  effects  (speckle)  in  COAT  systems. 

The  data  obtained  were  primarily  obtained  by  computer  simulations  using 
target  signatures  generated  by  GRC.  Comparisons  are  made  between 
analytical,  computer  simulation,  and  experimental  results. 

Section  IV  summarizes  the  principal  results  of  the  contract  and 
the  conclusions  that  can  be  drawn  from  the  work.  The  implications  for 
additional  research  that  should  follow  this  work  are  presented. 

An  appendix  contains  a development  of  the  theoretical  model  of 
the  effect  of  speckle  on  a COAT  system,  and  a statistical  model  of  a 
multidither  COAT  system.  In  combination,  these  analytical  models 
provide  a experimentally  verified  method  for  predicting  COAT 
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convergence  behavior  in  the  presence  of  speckle  modulation.  Computer 
simulation  results  using  artificial  speckle  modulation  are  also  in  good 
agreement  with  the  analytical  predictions.  This  appendix  forms  a basis 
for  understanding  the  computer  simulation  results  presented  in  Section  III. 


SECTION  II 


EXPERIMENTAL  TECHNICAL  ACCOMPLISHMENTS 


A.  EXPERIMENTAL  TARGET  SPECKLE  MEASUREMENTS 


In  the  first  part  of  the  experimental  study  reported  in  the  Interim 
Technical  Report,  we  characterized  the  general  properties  of  speckle 
patterns  associated  with  various  targets  to  identify  those  most  likely  to 
produce  degradation  in  COAT  system  performance.  Speckle  pattern 
measurements  were  obtained  on  targets  consisting  of  glass  spheres,  cones, 
and  cylinders  that  were  either  electroplated  with  nickel  or  painted  to  simulate 
the  effect  of  real  target  surfaces.  The  static  speckle  characteristics  were 

O 

measured  in  the  apparatus  shown  in  Figure  1.  The  output  of  a 6328  A helium- 
neon  laser  was  spatially  filtered  and  expanded  with  a standard  telescope. 

The  beam  was  truncated  at  its  full  width,  half  maximum  (FWHM)  by  an  iris, 
and  then  recollimated  to  0.  15  mm  diameter  (FWHM)  with  a second  telescope. 
The  beam  then  illuminated  the  target  at  a distance  of  about  two  meters  from 
the  exit  pupil.  Energy  backscattered  from  the  target  was  reflected  by  a mir- 
ror through  a pinhole  iris  into  a photomultiplier  (PMT)  detector.  Two  irises, 
0.  15  mm  and  0.  05  mm  in  diameter,  were  available  during  the  experiment  so 
that  speckle  data  could  be  recorded  when  the  receiver  and  transmitter  were 
of  equal  diameter  (D^  = = 0.  15  mm)  and  when  the  receiver  diameter 

was  equal  to  one  third  the  transmitter  diameter  (D^.  = 0.  15  mm, 

= 0.  05  mm).  The  target  stage  was  rotatable  and  the  PMT  was  trans- 
latable so  that  the  effects  of  target  rotation  and  translation  on  the  measured 
speckle  pattern  could  be  compared.  Standard  optical  chopper-synchronous 
detector  signal  processing  provided  .ood  signal-to-noise  (S/N)  ratios  even 
for  target  configurations  with  very  sn.  11  backscatter.  For  small  angle  tar- 
get rotation  or  small  receiver  translation,  the  speckle  patterns  obtained 
from  a given  target  were  identical.  Thus  continuous  target  rotation  was 
selected  for  the  remainder  of  the  measurements  made  during  the  program. 

A typical  speckle  pattern  for  metallized  sphere  is  shown  in 
Figure  2.  The  effect  of  aperture  averaging  of  the  speckle  data  can  be  seen 
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Figure  1„  Measurement  apparatus  for  study  of  static 
speckle  characteristics. 


by  comparing  the  upper  trace  taken  with  the  larger  aperture  to  the  lower 
trace  taken  with  the  smaller  aperture.  The  lower  trace  shows  higher  spatial 
frequency  content.  Similar  speckle  data  were  obtained  from  the  other 
metallized  targets  studied.  Maximum  local  contrast  ratios  as  high  as  0.  89 
and  as  low  as  0.  17  were  observed  under  various  experimental  conditions, 
with  most  of  the  data  showing  average  contrast  ratios  on  the  order  of  0.  5. 

Spatial  frequency  spectra  were  derived  from  the  speckle  data  of 
Figure  2 by  means  of  a digital  fast  Fourier  transform  (FFT)  routine.  These 
spectra  are  shown  in  Figure  3. 

In  the  Interim  Technical  Report  on  this  program  we  demonstrated 

1 2 

that  these  observed  spectra  agree  with  theory.  ’ The  dashed  curve  super- 
imposed on  each  figure  represents  the  calculated  effect  of  a one -dimensional 
aperture  integration  by  a slit  of  width  equal  to  the  aperture  diameter.  This 
procedure  permits  a qualitative  illustration  of  the  effects  of  aperture 
averaging.  As  expected,  the  spatial  frequency  spectrum  taken  with  the 
smaller  receiver  aperture  has  significantly  more  high  frequency  content. 

Rotation  of  the  targets  at  much  higher  speeds  than  those  used  in  the 
speckle  pattern  measurements  gave  real  time  speckle  modulations  which 
were  recorded  with  a spectrum  analyzer.  Figure  4(a)  shows  the  receiver 
noise  level  while  Figure  4(b)  shows  the  spectrum  of  the  speckle  modulations 
that  occur  with  a rotating  target.  The  FFT  data  in  Figure  4(b)  are  found 
from  the  spatial  frequency  data  of  Figure  3 by  using  the  relation 

f = zn^zf  , (1) 

t T s 

where  Z is  the  target  range,  f^.  is  the  temporal  frequency,  and  fg  is  the 
spatial  frequency.  As  can  be  seen,  the  spectrum  determined  by  the  two 
separate  techniques  are  in  good  agreement.  As  discussed  later,  except 
for  dither  modulations,  these  speckle  spectral  data  are  virtually  identical 
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Figure  3.  Spatial  frequency  spectra  of  receiver  modulations  generated  by 
rotation  of  a metallized  sphere  (same  as  Figure  2).  Dashed 
lines  show  theoretical  effect  of  aperture  averaging. 


Figure  4,  Speckle  spectra  for  smooth  spherical  metallized  target 
(a)  Stationary  target,  = 0.  (b)  Target  rotating  at  Sly  = 0.67  rps 

Spectrum  analyzer  bandwidth  is  10  Hz. 
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to  that  obtained  in  later  experiments  when  the  COAT  system  is  part  of  the 
optical  source. 

B.  COAT  SYSTEM  PERFORMANCE  IN  PRESENCE  OF  SPECKLE 
MODULATION 

In  general,  the  speckle  patterns  of  all  the  targets  studied  were  similar. 
Contrast  ratios  and  spectral  data  presented  in  the  previous  section  were  typi- 
cal for  metallic  targets.  These  targets  had  larger  contrast  ratios  than 
painted  targets  and  they  also  produced  more  intense  backs cattered  (180°) 
radiation.  Since  larger  contrast  ratio  speckle  will  produce  stronger  receiver 
modulations,  we  chose  a metallic  sphere  as  a representative  target  to  use 
with  the  COAT  system  to  study  interactive  effects.  The  stronger  backscatter 
also  helped  minimize  electronic  S/N  problems.  With  this  target,  we  observed 
the  convergence  behavior  as  a function  of  target  rotation  rate  and  of  the  ori- 
entation of  the  target  rotation  axis.  A diagram  of  the  experimental  apparatus 
is  shown  in  Figure  5. 

A description  of  the  COAT  system  portion  of  this  apparatus  is 
reported  in  reference  3.  The  output  beam  from  the  COAT  phasor  matrix, 
shown  expanded  in  the  figure,  passes  through  a reducing  telescope  and  a beam 
splitter.  The  primary  (stronger)  beam  then  passes  through  a turbulence  gen- 
erator and  a blooming  cell  placed  between  the  elements  of  a second  demagni- 
fying  telescope.  The  final  beam,  collimated  at  an  aperture  of  0.15  mm, 
illuminates  a rotatable  sphere  in  the  primary  target  plane.  The  backscattered 
energy  from  the  target  is  received  in  PMT  receiver  No,  1,  which  is  located 
close  to  the  target  demagnifying  lens. 

The  weaker  beam  from  the  beamsplitter  passes  through  an  identical 
demagnifying  telescope  before  illuminating  a fixed  secondary  target.  The 
backf  cattered  energy  from  this  target  is  received  by  a second  PMT,  located 
close  to  lens  f^.  The  two  optical  paths  are  thus  identical  except  for  the 
turbulence  and  blooming  distortions  and  the  tracking  and  focus  controls.  In 
the  interactive  experiments,  the  COAT  system  could  be  driven  from  either 
receiver  so  that  the  system  performance  with  and  without  speckle  effects 
could  be  compared,  An  iris  0.  15  mm  in  diameter  was  placed  in  front  of 
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receiver  No.  1 during  most  of  our  experiments  to  match  the  receiver  size  to 
the  diameter  of  the  transmitted  beam.  The  distance  from  the  target  to  both 
the  transmitter  and  receiver  was  one  meter. 

Since  it  was  not  feasible  to  place  a detector  on  the  rotating  target  to 
measure  the  convergence  level,  a second  beam  splitter  (not  shown)  was  used 
to  deflect  part  of  the  beam  illuminating  the  primary  target  into  a photodetector 
with  an  aperture  about  equal  to  the  null -null  converged  beam  diameter  at  the 
detector  location.  With  this  modification,  the  convergence  time  and  the  con- 
vergence level  of  the  beam  could  be  measured  while  the  target  rotated. 

We  first  performed  a series  of  convergence  ex.periments  with  the 
axis  of  rotation  of  the  target  oriented  at  6 = 0°,  45°,  and  90°  to  the  direction 
of  the  incident  laser  beam  to  simulate  several  typical  target/beam  encounters. 
The  receiver  spectrum,  the  average  converged  power,  and  the  COAT  system 
convergence  time  were  all  recorded.  The  blooming  and  turbulence  distortion 
generators  were  not  used  in  this  first  set  of  experiments. 

Figures  6 and  7 show  typical  experimental  data  for  the  two  cases  when 
the  angle  between  the  incident  beam  and  the  target  rotation  axis  is  90°  and  0°. 
The  45°  case  is  qualitatively  similar  to  the  90°  case  with  somewhat  less  degrad- 
ation observed.  In  Figure  6.  we  see  the  signals  recorded  at  three  target  rota- 
tion rates  for  the  90°  case.  The  upper  three  photographs  ere  the  receiver 
spectrum  and  the  lower  three  are  the  converged  power  on  the  target  (as  seen 
by  the  auxiliary  photodetector  discussed  above).  When  the  target  is  station- 
ary, the  normalized  average  power  is  defined  as  1 . 0 (no  degradation)  and 
the  system  convergence  time  is  2 msec.  The  receiver  spectrum  for  this 
case  is  simply  the  discrete  dither  spectrum. 

When  the  target  rotates  at  2tt  rad/sec,  the  speckle-induced  modula- 
tions cause  the  character  of  the  frequency  spectrum  to  change  from  a discrete 
to  a continuous  ■ pectrum  with  considerable  energy  falling  within  the  dither 
band.  The  average  converged  power  level  falls  to  0.47  and  the  time  history 
of  the  convergence  level  is  highly  oscillatory.  Rotation  of  the  target  at 
2tt  rad/sec  produced  a speckle  modulation  spectrum  which  covered  the  entire 
dither  band  and  produced  degradation  to  a power  level  of  0,  31.  In  these 
measurements  the  gain  was  kept  at  a high  level,  which  gave  optimum  con- 
vergence for  static  targets.  The  results  show  that  significant  convergence 
degradation  can  be  generated  by  speckle  modulations  from  moving  targets 
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Figure  6.  Experimental  convergence  data. 

Target:  Rough  metal  sphere,  90°  orientation  of  rotation  axis, 

^ - 0.488  pm. 

The  upper  pictures  are  outputs  of  a spectrum  analyzer  looking 
at  the  COAT  receiver  signal.  The  lower  pictures  show  the 
peak  target  irradiance  seen  by  a pinhole  detector  as  the  servo 
loop  is  closed. 

(a)  Target  rotation  = 0.  Average  power  after  convergence 
= 1. 

(b)  = tt  rad/sec  P = 0.70  (c)  firp  = 2-ir  rad/sec  P = 0.31. 
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Figure  7.  Effects  of  orientation  and  gain  in  speckle  experiments. 
Target  : Rough  metal  sphere,  X = 0.  488  pm,  target  rotation 
rate  £2-p  = 2t r rad/sec.  The  upper  pictures  are  outputs  of  a 
spectrum  analyzer  looking  at  the  COAT  receiver  signal.  The 
lower  pictures  show  the  peak  target  irradiance  seen  by  a pin- 
hole detector  as  the  servo  loop  is  closed,  (a)  Target  orienta- 
tion 9 = 0°,  static  optimization  (high  gain)  average  power  after 
convergence  = 0.92.  (b)  9 = 90°  static  optimization,  P = 0.31. 

(c)  9 = 90°  dynamic  optimization  (low  gain)  P = 0. 93. 


when  the  COAT  system  gain  is  optimized  for  a static  target.  If  the  gain  is 
reduced,  the  degradation  is  less  severe  because  the  servo  bandwidth  in  each 
channel  is  reduced.  Thus  the  overall  convergence  level  increases  because 
less  speckle  modulation  is  passed  by  the  servo  system  synchronous  detectors. 
We  show  these  effects  in  Figures  7(b)  and  7(c),  with  the  axis  of  rotation  at 
90°,  and  a rotation  rate  of  2ir  rad/sec.  Decreasing  the  static  optimization 
gain  level  by  -15  dB  produces  a dynamic  convergence  optimization  in  which 
the  mean  convergence  level  is  0.  93,  compared  with  0.  31  for  the  high-gain 
case. 

In  Figure  7(a),  with  the  target  rotation  axis  oriented  parallel  to  the 
beam  (i.  e.  , 0°),  we  observe  very  little  convergence  degradation.  Here  we 

see  for  high  gain  (static  optimization)  only  minimal  convergence  degradation 
to  a level  of  0,  92  at  a target  rotation  speed  of  2tt  rad / sec. 

Figures  8 and  9 summarize  the  observed  convergence  levels  and 
convergence  times  for  all  the  interactive  experiments.  Referring  to  the  0° 
orientation  case  in  Figure  8(a),  the  system  performance  was  first  optimized 
with  the  target  stationary;  the  convergence  level  was  then  recorded  as  a func- 
tion of  rotation  rate.  For  this  case,  optimization  of  the  system  with  the  tar- 
get rotating  produced  the  same  results.  The  small  variations  in  converged 
power  level  for  different  rotation  rates  are  due  to  glint  tracking  by  the  COAT 
system,  which  caused  the  converged  beam  to  partially  miss  the  detector 
aperture.  Visual  observation  of  the  converged  beam  profile  on  a television 
monitor  showed  almost  no  degradation. 

With  the  rotation  axis  at  90  degrees  (Figure  8(b)),  the  largest 
degradation  occurs  with  high  loop  gain  (static  optimization)  and  the  smallest 
degradation  occurs  with  low  loop  gain  (dynamic  optimization).  The  inter- 
mediate or  45°  orientation  case  in  Figure  8(c)  shows  speckle  convergence 
degradation  to  a level  of  0.  45  with  static  optimization  and  recovery  to  a level 
of  0.  80  with  dynamic  optimization.  The  average  convergence  level  for  all 
three  orientations  under  static  and  dynamic  optimization  conditions  is  shown 
in  Figure  8(d). 


Minimum  convergence  time  and  maximum  convergence  target  irradiance 
achieved  by  adjusting  servo  loop  gain. 
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Figure  8.  Summary  of  observed  convergence  levels  in  the 
presence  of  speckle  effects,  (a)  Angle  between 
t rget  rotation  axis  and  incident  beam,  9 = 0°. 

( ) 3 = 90°.  (c)  0 = 45°.  (d)  Data  of  (a),  (b), 

and  (c)  averaged. 
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Figure  9.  Summary  of  convergence  times  in  the  presence  of 
speckle  modulations  (0  = orientation  of  rotation 
axis  relative  to  incident  laser  beam). 

(a)  0 = 0°,  (b)  0 = 90°,  (c)  0 = 45°,  (d)  average 
of  (a),  (b),  and  (c). 
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In  addition  to  improving  the  overall  convergence  level,  another  effect 
of  reducing  the  COAT  system  loop  gain  is  to  increase  convergence  time. 

The  observed  convergence  times  are  shown  in  Figure  9.  Typical  convergence 
times  for  static  targets  are  1 to  2 msec.  The  convergence  time  for  0°  target 
orientation  is  essentially  the  same  at  all  rotation  speeds  because  there  was 
little  convergence  degradation,  and  in  all  cases  since  the  loop  gain  is  set  to 
its  maximum  value.  With  dynamic  optimization  of  the  convergence  level  (or 
lower  loop  gain),  the  convergence  time  increases  substantially,  as  shown  in 
Figures  9(b)  and  9(c).  For  these  measurements,  convergence  time  is  defined 
as  the  time  for  the  target  irradiance  to  rise  from  10%  to  90%  of  the  final 
average  level.  For  the  "rotation  optimization"  data,  the  gain  was  adjusted 
for  maximum  convergence  at  each  rotation  velocity,  resulting  in  a different 
convergence  time  for  each  measurement.  This  effect  is  not  a direct  con- 
sequence of  the  target  rotation  velocity,  but  rather  of  the  system  gain 
required  to  produce  dynamic  optimization  at  each  velocity. 

C.  COAT  SYSTEM  PERFORMANCE  IN  THE  PRESENCE  OF  SPECKLE 
MODULATION,  BLOOMING,  AND  ARTIFICIAL  TURBULENCE 

Experimental  thermal  blooming  and  turbulence  compensation  have 
been  studied  with  speckle  modulation  effects  also  present.  The  apparatus 
used  in  these  studies  is  identical  to  that  used  in  earlier  experiments  which 
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had  no  speckle  effects.  Measurements  were  performed  using  a flowing 

liquid  blooming  cell,  with  the  cell  occupying  the  first  81%  of  the  focused, 

22  cm  propagation  path.  This  condition  provides  better  correction  than  the 

4 

case  where  the  cell  is  in  the  last  81%  of  the  propagation  path. 

To  ensure  that  the  blooming  distortions  as  a function  of  propagation 
distance  are  the  same  in  the  laboratory  at  0.488  [im,  the  blooming  strength, 
the  total  absorption,  a Z,  and  the  Fresnel  number  must  be  held  constant.  The 
absorption  coefficient  is  a and  the  total  propagation  distance  in  the  blooming 
medium  is  Z. 


The  F resnel  number,  is  conveniently  defined  by 


where  k - Ztt /\  and  \ is  the  wavelength,  and  D^,  is  the  transmitter  diameter. 
The  blooming  strength  is  measured  by  the  product  oT  . P^,  where 

J_j  i 

P„  is  the  laser  power  at  the  entrance  to  the  blooming  medium  and  o'  is 
defined  for  a gaseous  medium  by 
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The  quantities  yq>  Y>  and  V are  the  molecular  polarizability  of  the  gas, 
the  specific  heat  (y  = Cp/C^.),  the  static  gas  pressure,  and  the  transverse 
wind  velocity.  For  a liquid,  o is 

i-J 


(%).. 
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where  the  liquid  density  is  p,  the  nominal  index  of  refraction  is  n^,  and 
dn/dT  is  the  change  of  index  with  temperature. 

Because  the  net  blooming  cell  transmission  (windows  plus  liquid)  was 
only  0.  1,  the  receiver  aperture  had  to  be  chosen  with  a diameter  twice  the 
transmitter  in  order  to  achieve  an  adequate  S/N  level.  Some  aperture 
averaging  over  the  data  in  Section  II-B  is  thus  occurring.  With  the  speckle- 
generating  target  in  motion,  system  convergence  was  unstable  at  all  but  the 
highest  power  levels  with  a smaller  receiver,  the  instability  being  produced 
more  by  the  low  S/N  than  by  speckle  effects.  By  increasing  the  receiver 
aperture  to  twice  the  transmitter  size,  we  were  able  to  obtain  more  reliable 
and  reproducible  data  as  discussed  in  the  following  subsections. 
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1. 


Blooming  and  Speckle 


In  this  series  of  experiments,  scaled  according  to  Eqs.  (2)  and  (3), 
we  wanted  to  observe  any  possible  interactions  between  the  speckle  and 
blooming  phenomena  and  how  such  interactions  might  affect  COAT  perfor- 
mance. The  blooming  cell  shown  schematically  in  Figure  5 was  filled  with 
a methanol  and  iodine  solution  which  gave  a transmission  factor  of 
exp  [-ah)  = 0.3  for  the  liquid  alone  [ah  = 1.  2).  The  cell  was  operated  at 
a transverse  wind  velocity  between  0.  15  cm/ sec  and  0.  55  cm/sec  which 
produced  significant  convection-dominated  thermal  blooming.  The  trans- 
mitter diameter  was  0.  1 cm.  The  target  rotation  speed  selected  for  the 
speckle  measurements  was  always  2ir  rad/sec,  the  speed  at  which  maximum 
convergence  degradation  is  observed  with  the  COAT  system  in  the  high  loop 
gain  condition  (static  optimization,  see  previous  section).  Observations  of 
the  convergence  level  were  made  at  various  laser  powers  with  and  without 
speckle  effects,  with  and  without  active  COAT  phase  compensation;  and  at 
both  the  high  and  low  servo  loop  gain  settings.  The  results  of  these  experi- 
ments are  shown  in  Figure  10. 

In  the  absence  of  speckle,  the  COAT  system  improves  the  peak  target 
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irradiance  by  a factor  of  2.  3 and  increases  the  critical  transmitter  power 
by  roughly  the  same  factor.  The  improvements  in  the  target  irradiance  at 
the  corrected  critical  power  (P^  « 3.  5)  is  roughlv  a factor  of  3.  1. 

In  the  presence  of  speckle  effects,  with  the  system  optimized  for 
best  correction  in  the  absence  of  speckle  (high  gain),  rotation  of  the  target 
at  2tt  rad/sec  produces  a converged  irradiance  degradation  of  0.  52  of  the 
maximum  at  the  critical  power.  Reduction  of  the  servo  loop  gain  by  about 
10  to  15  dB  raises  the  converged  irradiance  back  to  a level  of  0.  81  of  maxi- 
mum. This  is  qualitatively  similar  to  what  we  observed  in  the  nonblooming 
experiments.  Without  the  blooming,  however,  we  observed  a more  severe 
speckle  degradation  and  a better  recovery,  from  0.  31  to  0.  93  of  the  free- 
space,  no-speckle  maximum. 


Transmitter  power  at  which  the  target  irradiance  is  maximum  when 
blooming  effects  are  present. 
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PT  = TRANSMITTER  POWER,  ARBITRARY  UNITS 

Figure  10.  Target  (focal-plane)  irradiance  with  blooming 

and  speckle,  with  and  without  COAT  correction. 
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Nevertheless,  at  high  power  levels,  the  converged  irradiance  is 
higher  with  COAT  correction  than  without  in  all  cases,  It  is  interesting  to 
note,  however,  that  at  lower  power  levels  (P^,  < 2.  0 in  Figure  10),  the 
speckle  effects  dominate  the  blooming  so  that  the  high  gain,  COAT-corrected 
case  actually  produces  a slightly  lower  target  irradiance  than  the  no-COAT 
case. 


2.  Blooming,  Artificial  Turbulence,  and  Speckle 

4 

On  a previous  contract,  we  established  experimentally  that  turbulence 

distortions  could  be  removed  by  a COAT  system,  even  though  the  absolute 

Strehl  ratio  was  degraded  by  blooming.  Thus  the  nonlinear  medium  did  not 

limit  the  COAT  system's  ability  to  correct  for  turbulence. 

In  this  program,  the  ability  of  the  COAT  system  to  simultaneously 

compensate  turbulence  and  blooming  was  also  studied,  but  with  the  addition 

of  moving-target  speckle  noise  as  a new  parameter.  Seven  cases  were 

studied,  with  blooming  present  in  each  one,  as  shown  in  Table  2.  By  activat- 
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ing  the  artificial  turbulence  generator  (described  and  evaluated  in  the  pre- 

4 

vious  contract  ),  we  were  able  to  study  the  combined  effects  of  blooming  and 
turbulence  on  COAT  system  convergence  in  the  presence  and  absence  of 
speckle  modulation.  The  results  without  speckle  modulation  are  shown  in 
Figure  11  and  with  speckle  modulation  in  Figure  12. 


PT  = TRANSMITTED  POWER,  ARBITRARY  UNITS 


• ••  . s. 

Figure  11.  COAT  convergence  with  blooming  and  turbulence, 
but  with  no  speckle  modulation  effects  (stationary 
target).  The  Cpj  = 0 curves  are  the  same  as  in 
Figure  10. 
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PEAK  TARGET  IRRADIANCE.  ARBITRARY  UNITS 


Figure  12.  COAT  convergence  with  blooming,  artificial 
turbulence,  and  speckle  modulation  effects. 
Target  rotation  rate,  = 6.3  rad/sec. 
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There  are  several  things  to  note  about  the  data  in  Figure  li.  First, 

the  COAT  system  cannot  fully  correct  for  the  turbulence;  the  strength  of  the 

turbulence  is  too  great  for  the  finite  number  of  phase  correction  elements. 

The  effect  is  seen  at  low  transmitter  power  where  the  slope  of  the  "COAT 

corrected,  C..  t 0"  curve  falls  below  the  free  space  line  and  the  "NO  COAT, 
2 ^ 

CN  = 0"  line.  Significant  turbulence  compensation  is  occurring,  however, 

since  the  "COAT  corrected,  CN^  t 0"  curve  falls  well  above  the  "NO  COAT, 

C^  t 0"  curve.  Second,  in  the  presence  of  blooming  and  turbulence,  the 

COAT  system  produces  about  a factor  of  2.  5 increase  in  peak  target  irra- 

diance.  In  fact,  at  all  power  levels  more  than  a factor  of  2 correction  is 

4 I 

observed.  Finally,  we  conclude,  as  we  have  previously,  that  the  turbulence 
compensation  capacity  of  a multidither  COAT  system  is  not  significantly 
affected  by  the  presence  of  blooming  distortions. 

When  speckle  effects  are  added  to  these  experiments,  a reduction  in 
peak  target  irradiance  occurs  for  the  COAT -correr m f cases,  as  expected. 
The  results  are  shown  in  Figure  12.  All  of  the  COAT-corrected  data  in 
Figure  12  were  taken  with  a lowered  servo  loop  gain  ("dynamic  optimization" 
as  used  earlier).  Very  little  turbulence  compensation  is  observed;  the  con- 
vergence degradation  caused  by  speckle,  about  a factor  of  2,  is  just  enough 
to  offset  the  improvement  in  Figure  11  resulting  from  turbulence  compensa- 
tion, also  about  a factor  of  2. 

Note,  however,  that  blooming  compensation  is  only  slightly  affected 

2 

by  speckle  as  seen  earlier  in  Figure  10.  With  no  turbulence  (C = 0),  the 
peak  irradiance  is  increased  by  COAT  correction  about  a factor  of  1.  7,  which 
is  less  than  the  factor  of  2.  3 in  Figure  11,  but  still  significant.  The  lack  of 
turbulence  compensation  observed  with  these  severe  speckle  effects  is  not 
well  understood  since  significant  blooming  compensation  occurs,  indicating 
that  the  COAT  system  phase  correction  capabilities  are  not  completely  lost. 
As  a possible  explanation,  we  hypothesize  that  the  speckle  degrades  the  high 
frequency  error  correction  capability  of  the  COAT  servo  more  than  the  low- 
frequency  correction.  In  effect,  the  effective  servo  error  correction  band- 
width has  been  reduced.  The  very  low  electronic  S/N  ratios  present  in  our 
experiments  may  also  contribute  to  the  poor  turbulence  compensation  with 
speckle.  That  is,  when  the  S/N  is  low,  any  noise-like  disturbances  (such  as 
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speckle)  have  more  effect  on  the  system  and  the  high  frequency  servo  error 
correction  capability  is  reduced.  We  have  observed  similar  poor  turbulence 
correction  behavior  in  other  experiments  where  there  were  no  speckle  effects 
but  which  had  very  poor  electronic  S/N.  At  this  writing,  however,  we  have 
no  way  to  uniquely  isolate  electronic  S/N  effects  from  purely  speckle-induced 
effects  for  the  data  in  Figure  12.  The  data  definitely  indicate  that  the 
turbulence  compensation  capability  of  a COAT  system  i3  strongly  affected  by 
speckle,  and  that  the  blooming  correction  is  also  reduced,  but  much  less 

so. 
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SECTION  III 


ANALYTICAL  AND  COMPUTER  SIMULATION 
TECHNICAL  ACCOMPLISHMENTS 


A.  THEORETICAL  ANALYSIS  OF  MULTIDITHER  COAT  SERVO 

RESPONSE  TO  SPECKLE  MODULATION 

A comprehensive  analysis  has  been  developed  during  the  course  of 
this  program  to  explain  the  convergence  behavior  of  COAT  systems  in  the 
presence  of  spurious  modulations  due  to  speckle.  These  results  have  been 
partially  reported  in  Interim  Report  No.  1 of  this  contract.  An  updated 
documentation  of  these  results  is  given  in  the  Appendix  of  this  report  and 
summarized  here  very  briefly. 

The  analysis  predicts  average  convergence  levels  as  a function  of  the 
average  power  in  the  spurious  electrical  signal  in  the  COAT  servo  due  to 
speckle-induced  intensity  modulations  at  the  receiver.  In  order  to  predict 
the  convergiiiice  level,  we  must  assume  that  the  average  power  spectrum  of 
the  modulations  is  known  or  can  be  determined.  Since  the  power  spectrum  is 
a function  of  the  degree  of  convergence  of  the  incident  beam,  of  target  sur- 
face characteristics,  and  of  target  motion,  all  this  information  must  be  known. 
The  analysis  neglects  the  change  in  power  spectrum  with  convergence  level. 
Nevertheless,  in  order  to  be  conservative  we  have  assumed  a target  illum- 
ination near  full  convergence  since  this  will  give  speckle  lobes  on  the  order 
of  a transmitter  diameter.  Lower  convergence  levels  will  result  in  smaller 
speckles  which  can  be  more  easily  averaged  by  the  receiver  aperture.  We 
are  thus  considering  only  power  spectrums  which  are  worst  cases  in  virtually 
all  physical  applications.  Knowing  the  power  spectrum  and  the  COAT  servo 
bandwidth  allows  us  to  calculate  the  average  power  in -the  spurious  electronic 
signals.  Using  a statistical  model  of  the  convergence  process,  we  then  cal- 
culate the  average  power  in  the  true  error  signal  in  the  COAT  servo  at  the 
various  dither  frequencies.  The  average  convergence  level  is  then  determined 
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based  on  the  assumption  that  the  power  ratio,  p\  of  the  dither  signal  to  the 
speckle -induced  signal  after  demodulation  is  unity.  Computer  simulation 
studies  discussed  in  the  next  section  have  validated  ihis  assumption. 

In  a typical  multidither  COAT  system,  the  servo  bandwidth  in  a given 
channel  is  a function  of  the  gain  in  that  channel.  The  gain,  on  the  other  hand, 
depends  on  the  overall  convergence  level  and  the  error  in  that  channel.  The 
average  power  in  the  speckle  signal  will  therefore  also  vary  as  a function  of 
convergence  level  due  to  this  variation  in  bandwidth.  This  effect  can  be 
included  in  the  analysis  by  calculating  an  average  bandwidth  for  all  the  chan- 
nels as  a function  of  convergence  level.  This  would  involve  another  inde- 
pendent analysis  of  the  COAT  servo  and  because  of  limited  time  was  not  done 
as  part  of  this  study.  In  lieu  of  a complete  first-principles  calculation,  we 
have  verified  the  analysis  by  comparing  its  predictions  to  computer  simula- 
tions using  an  artificial  modulation  function  whose  power  spectrum  consisted 
of  very  narrow  bands  around  each  dither  frequency.  The  bands  were  chosen 
narrow  enough  so  that  they  would  always  be  narrower  than  the  servo  band- 
width. Thus,  there  could  be  no  variation  in  the  spurious  modulation  power 
due  to  changes  in  the  convergence  level.  The  results,  which  show  excellent 
agreement  between  analysis  and  computer  simulation,  are  shown  in  Figure  7 
of  the  Appendix.  If  the  analysis  of  the  bandwidth  variations  is  completed  as 
part  of  future  programs,  we  expect  to  be  able  to  show  similar  agreement 
between  analysis  and  computer  simulation  for  any  realistic  power  spectrum. 

The  principal  result  of  this  analysis  is  that  we  can  now  analytically 
predict  the  average  convergence  level  of  a COAT  system  if  two  parameters 
are  well  specified:  power  spectrum  of  the  speckle  modulations  and  COAT 

servo  bandwidth.  If  these  variables  can  be  specified  as  a function  of  con- 
vergence level,  then  iterative  calculations  can  be  made  which  will  predict 
the  correct  level.  If  they  cannot  be  specified  as  a function  of  convergence 
level,  then  less  accurate  but  conservative  predictions  can  still  be  made. 
Another  important  result  of  the  analysis  is  that  comparison  with  computer 
simulations  indicate  that  the  maximum  fluctuations  about  the  average  con- 
vergence level  are  also  predictable.  The  computer  simulations  (as  well  as 
the  experimental  observations  discussed  in  the  previous  chapter)  showed  that 
the  convergence  level  in  the  presence  of  speckle  effects  oscillates  randomly 
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between  two  levels.  The  analysis  indicates  that  the  maximum  level  is 
reached  when  the  power  in  the  speckle  signal  is  twice  the  power  in  the  dither, 
that  is  when  p1  = 0,  5,  and  the  minimum  level  occurs  when  p'  = 2.  0. 

B.  COMPUTER  SIMULATION  RESULTS 

During  the  first  phase  of  this  contract,  computer  simulations  were 
run  using  a set  of  three  realistic  speckle  modulation  data  generated  by  the 
General  Research  Corporation  (GRC).  ^ The  physical  scenario  consisted  of 
a rotating  target  sphere  (1  m radius)  uniformly  illuminated  with  10.  6 pm 
coherent  radiation  over  a 10  cm  diameter  spot  at  a range  of  2 km  from  the 
COAT  receiver.  This  illumination  function  corresponds  roughly  to  a 
diffraction-limited  spot  from  a 25  cm  diameter  transmitter.  The  data  itself 
consisted  of  time  histories  of  the  intensity  variations  at  the  receiver  for 
three  target  rotation  rates,  D=  0,01,  2,  and  10  rad/sec,  with  the  rotation 
axis  oriented  perpendicular  to  the  illuminating  beam  direction.  These  results 
along  with  some  of  our  reservations  about  the  data  were  reported  in  the  last 
interim  report. 

Since  the  previous  report,  we  have  resolved  our  questions  about  their 
model  and  a new  set  of  data  was  supplied  by  GRC.  This  data  included  aper- 
ture averaging  over  an  annular  receiver  surrounding  a transmitter  with  an 
outside  diameter  of  1.2m  and  an  inner  diameter  of  1,0m,  A new  calculation 
technique  using  statistical  methods  to  calculate  the  intensity  variations  or 
modulation  functions  was  also  employed.  This  technique  overcomes  the 
problem  encountered  in  the  earlier  data  where  the  modulations  were  periodic 
with  a period  of  approximately  0,  5 msec  caused  by  a finite  number  of  sample 
points  in  a point  target  model.  The  new  data  again  consisted  of  three  modu- 
lation functions  representing  the  same  physical  scenarios,  except  that  the 
slowest  target  rotation  rate  was  = 0.4  rather  than  0.  01  rad/sec,  A 
schematic  diagram  of  the  scenario  is  shown  in  Figure  13. 

We  found  that  all  three  of  these  modulation  functions  produced 
noticeable  degradation  in  COAT  performance.  The  amount  of  degradation 
depends  on  both  the  speckle  modulation  itself  and  on  the  choice  of  servo 
parameters.  The  principal  observable  effects  are  (1)  a decrease  in  the 
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mean  convergence  level,  and  (2)  random  oscillation  of  the  mean  level  at 
frequencies  on  the  order  of  0.  5 to  1.0  kHz.  Of  the  various  servo  variables, 
the  open  loop  gain  and  the  signal  amplitude  limiter  (clipper,  see  Fig- 
ure 14)  had  the  most  noticeable  influence  in  reducing  these  effects.  Both 
the  experiments  and  the  computer  simulations  demonstrate  that  reducing 
from  its  optimum  value  for  performance  without  speckle  will  increase  the 
degraded  mean  convergence  level  when  speckle  modulations  are  present. 

This  has  the  undesirable  effect  of  also  increasing  the  convergence  time 
(decreasing  the  servo  control  bandwidth).  Signal  amplitude  limiting  can  also 
produce  an  increase  in  the  mean  convergence  level  for  selected  values  of  the 
amplitude  limit  and  for  certain  ranges  in  the  value  of  G^.  The  effect  on 
convergence  time,  on  the  other  hand,  is  much  less  noticeable. 

In  order  to  pick  the  best  values  of  these  parameters,  an  optimization 
study  was  conducted  using  the  realistic  speckle  modulation  function  with  the 
target  rotating  at  2.  0 rad/sec,  chosen  because  it  produced  the  most  severe 
degradation  in  COAT  performance.  Figures  15  and  16  summarize  the  data 
gathered  during  this  study.  Figure  15  shows  that  for  each  clipping  level  there 
is  an  optimum  gain  where  the  convergence  time  is  minimum.  Increasing  or 
decreasing  the  gain  from  this  level  causes  slower  convergence.  In  Figure  16, 
for  the  same  set  of  clipping  levels,  we  see  that  the  highest  mean  convergence 
levels  occur  at  lower  gains  than  those  associated  with  the  smallest  conver- 
gence times.  The  bars  in  this  figure  represent  the  extremes  of  the  oscillatory 
variation  in  the  convergence  level  referred  to  earlier.  The  value  of  signal 
clipping  is  demonstrated  by  comparing  Figure  16(a)  and  16(d).  In  Figure  16(d) 
the  low  clipping  greatly  reduces  the  oscillations  in  the  convergence  level  and 
improves  the  mean  convergence  level,  except  at  very  low  gains.  In  Figure 
15(d),  however,  we  see  that  this  value  of  clipping  level  tends  to  increase 
the  convergence  time.  In  the  optimization  study,  these  tradeoffs  were  con- 
sidered, and  we  chose  G^  = 1000  and  a clipping  level  of  ±0.3  as  a reasonable 
compromise  for  both  high  convergence  level  and  fast  convergence  time.  All 
three  of  the  GRC  scenarios  were  then  simulated  using  these  values  of  clipping 
and  gain. 
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Figure  16.  Summary  of  convergence  level  optimization  data 
in  COAT  simulations. 


Figures  17,  18,  and  19  show  the  three  GRC  modulation  functions,  their 
Fourier  spectrum,  and  the  resulting  COAT  performance  in  the  form  of  a time 
history  of  the  intensity  on  a target  glint  normalized  to  unity  at  full  con- 
vergence. (The  two  time  histories  superimposed  on  each  other  correspond 
to  COAT  performance  with  and  without  speckle  modulations  for  easy 
comparison.  ) 

The  most  obvious  conclusion  from  these  simulations  is  that  the 
modulation  power  in  the  dither  band  is  the  overriding  speckle  parameter 
which  will  affect  COAT  performance.  At  a slow  rotation  rate  (Figure  17),  most 
of  the  power  is  at  the  frequencies  below  the  dither  band.  At  the  higher  rota- 
tion rate  of  2 rad/sec  (Figure  18),  much  of  the  power  spectrum  falls  in  the 
dither  band.  At  the  highest  rotation  rate,  10  rad/sec  (Figure  19),  the  power 
spectrum  spreads  over  a large  frequency  range;  consequently  there  is  less 
power  in  the  dither  band  than  in  Figure  19.  The  most  severe  degradation 
occurs  for  the  case  in  Figure  18  where  most  of  the  spectral  power  is  in  the 
dither  band.  Increasing  or  decreasing  the  rotation  rate  redistributes  the 
power  spectrum  outside  the  dither  band  with  a resultant  improvement  in  COAT 
performance . 

These  results  can  be  compared  to  analytical  predictions  subject  to 
the  limitations  described  earlier.  At  this  writing,  we  can  only  estimate  what 
the  average  servo  bandwidth  corresponding  to  each  simulation  actually  is, 
since  the  actual  bandwidth  is  a function  of  time  and  varies  from  channel  to 
channel.  Our  best  estimate  indicates  that  an  average  of  ±100  Hz  around  each 
dither  frequency  is  a reasonable  number  for  all  three  simulations,  for  the 
choice  of  Gl  = 1000  and  clipping  level  = ±0.3.  We  thus  calculate  the  three 
speckle  coefficients*  based  on  this  bandwidth.  The  v<  ults  of  these  calcula- 
tions along  with  the  corresponding  predicted  convergence  levels  using  Fig- 
ure A-9  of  the  appendix  are  shown  in  Table  3.  The  observed  convergence 
levels  in  the  computer  simulations  are  also  shown  for  comparison.  Note  that 
the  agreement  is  reasonable,  with  the  analysis  tending  to  be  pessimistic  in  its 
prediction.  With  a proper  calculation  of  an  average  servo  bandwidth,  we 
feel  the  agreement  will  improve. 


’'See  appendix  for  definition  of  speckle  coefficient,  Cg. 

49 


5045-11 


1.0 


uj  0.8 
> 

LU 

-J 

UJ  0.6 
O 

<c>  2 

g 0.4 

UJ 

> 

O 0-2 
a 

0 


Figure  17.  Summary  of  computer  simulation  results 
with  realistic  speckle  modulations  from  a spherical 
target  rotating  at  0.4  rad/sec  at  a range  of  2 km  from 
receiver,  (a)  Fourier  spectrum  of  speckle  modulations 
(b)  Time  history  of  speckle  modulations,  (c)  Time 
history  of  GOAT  performance. 
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Figure  18,  Summary  of  computer  simulation 
results  with  realistic  speckle  modulations 
from  a spherical  target  rotating  at  2.0  rad/sec 
at  a range  of  2 km  from  receiver,  (a)  Fourier 
spectrum  of  speckle  modulations,  (b)  Time 
history  of  speckle  modulations,  (c)  Time 
history  of  COAT  performance. 


''  .•r<.vT,'  i',j 


>;sm 


0 8 


5046-16 


uj  0 6 

Q 

=> 

(a)  □ 0,4 

CL 

2 

< 0.2 


k H DITHER  BAND 


0 20  40  60  80  1 00  1 20 

FREQUENCY,  khz 


5045-17 


5045-18 


Figure  19.  Summary  of  computer  simulation 
results  with  realistic  speckle  modulations  from 
a spherical  target  rotating  at  10.0  rad/ sec  at  a 
range  of  2 km  from  receiver,  (a)  Fourier  spec- 
trum of  speckle  modulations  , (b)  Time  history 

of  speckle  modulations,  (c)  Time  history  of 
COAT  performance. 
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TABLE  3.  COMPARISON  OF  ANALYTICAL  SPECKLE  MODEL 
WITH  COMPUTER  SIMULATIONS  FOR  REALISTIC  SPECKLE 

MODULATIONS 


Target 

Rotation 

Rate 

(rad/ sec) 

Computer  Simulation 
Convergence  Level 

Calculated 
fr  om 

Simulation 

Analytically  Predicted 
Convergence  Level 

max 

min 

avg 

max 

min 

avg 

0 

0.  95 

0.  95 

0.  95 

0 

0.  95 

0.  95 

0.  95 

0.  4 

0.  86 

0.  72 

0.  74 

0.  07 

0.  87 

0.  68 

0.  79 

2 

0.  42 

0.  22 

0.  32 

0.  24 

0.  39 

0.  05 

0.  18 

10 

0.  62 

0.  40 

0.  51 

0.  18 

0.  56 

0.  15 

0.  36 

Other  servo  parameters  can  also  be  used  to  some  degree  to  help 
improve  COAT  performance.  Increasing  the  dither  amplitude,  , for  example, 
increases  the  dither  signal  power  relative  to  the  speckle  modulation  power, 
with  a resulting  increase  in  convergence  level.  Beyond  a given  amplitude 
which  depends  on  the  severity  of  the  degradation,  however,  one  reaches  a 
point  of  diminishing  return  because  the  maximum  possible  convergence  level 
decreases  with  increasing  dither  amplitude.  This  result  is  illustrated  in 
Figure  20.  The  solid  curve  is  the  theoretical  maximum  convergence  level, 
while  the  data  points  are  from  computer  simulation  runs  using  the  artificial 
speckle  modulation  function  of  Eq.  A- 27,  in  the  appendix,  with  the  amplitudes 
a.,  chosen  randomly  between  0 and  0.  012.  This  corresponds  to  a speckle 
coefficient,  Cg  = 0.095.  If-  4*  = 20°  is  considered  the  nominal  dither  amplitude, 
one  can  achieve  some  improvement  in  convergence  level  by  increasing  ^ to 
35°,  but  beyond  this  value  increasing  ^ will  only  cause  more  degradation. 

During  this  program  we  also  found  that  changing  to  square  wave 
demodulation  in  the  synchronous  detection  process  seemed  to  reduce  the 
convergence  level  degradation  caused  by  speckle  modulations.  The  observed 
effect  is  shown  in  Figure  21.  The  two  time  histories  have  identical  servo 
parameters  except  that  part  a)  has  square  wave  demodulation  while  part  b) 
has  sine  wave  demodulation.  The  limited  time  available  for  the  study  of 
square  wave  demodulation  only  permitted  observation  of  this  effect  in  a few 
simulations.  We  therefore  offer  no  explanation  as  to  why  it  should  produce 
this  apparent  improvement.  We  only  document  it  in  this  report  for 
completeness. 
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The  principal  accomplishments  of  the  analytical  study  are  twofold. 
First,  it  has  provided  insight  and  a good  understanding  of  how  speckle- 
induced  modulations  can  interfej  with  normal  COAT  servo  operations  and 
how  the  servo  itself  responds  eir  presence.  Second,  we  have  developed 

a tool  which  predicts  the  convergence  level  of  a COAT  array  subject  to 
speckle  degradation,  given  that  the  power  spectrum  of  the  speckle  modula- 
tions and  the  servo  bandwidth  is  known.  The  predictions  correspond  very 
well  both  qualitatively  and  quantitatively  to  the  results  observed  in  the  com- 
puter simulations.  However,  it  is  difficult  to  adequately  compare  the 
theoretical  analysis  to  the  experimental  data,  because  of  the  difficulty  in 
measuring  the  synchronously  detected  speckle  modulation  power  in  the  COAT 
servo.  If  careful  measurements  are  taken  and  if  the  variation  of  servo  band- 
width as  a function  of  convergence  level  is  included  in  the  analysis,  we  expect 
the  analytical  predictions  to  be  just  as  good  for  experimental  data  as  they  are 
for  simulation  data. 

C.  COMPARISON  OF  EXPERIMENTAL  AND  COMPUTER  SIMULATION 

RESULTS 

As  noted  previously,  the  computer  simulation  results  and  the  experi- 
mental data  agree  quite  well  qualitatively.  For  quantitative  comparison,  it 
is  necessary  to  scale  the  computer  simulation  parameters  to  the  laboratory 
experiment.  Unfortunately,  the  simulation  and  experimental  parameters  were 
chosen  independently  and  consequently,  the  principal  cases  studied  do  not 
correspond  exactly.  In  particular,  effects  such  as  aperture  integration  and 
target  illumination  function  were  not  included  in  the  scaling  and  these  effects 
can  significantly  alter  the  spatial  and  temporal  frequency  spectra  generated 
by  target  speckle  motion.  Consequently,  differences  in  speckle  degradation 
of  convergence  will  be  observed. 

In  spite  of  these  limitations,  it  is  still  useful  to  compare  two  cases 
whose  Fourier  power  spectrum  are  very  similar.  Throughout  this  report  we 
have  stressed  that  speckle-induced  power  in  the  dither  band  is  the  crucial 
parameter  in  determining  COAT  performance.  We  will  now  demonstrate 
that  computer  simulations  and  experimental  data  both  agree  quantitatively 
as  well  as  qualitatively  when  we  compare  cases  subjected  to  the  same  speckle 
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power  spectrum.  For  a rigorous  comparison  of  the  rotating  sphere 
scenarios  used  in  the  experiments  and  simulations,  one  should  first  scale 
the  diameter  of  transmitter,  receiver,  and  illuminated  target  spot  size  to 
have  the  same  ratios;  second,  we  should  require  the  contrast  ratios  to  be 
the  same;  and  third,  the  spatial  frequency  distributions  should  be  proportional. 
With  all  the  above  requirements  met,  we  would  then  scale  the  target  range  and 
rotation  rate  using  the  relationship 


n z 

T 


constant  , 


where  is  the  target  rotation  rate,  Z is  the  target  range,  and  is  the 
transmitter  diameter.  This  insures  that  the  temporal  frequency  distribution 
is  the  same.  Since  the  analytical  model  used  at  GRC  for  the  speckle  modula- 
tion calculations  assumed  random  spatial  phase  distribution  for  the  reflected 
target  radiation,  and  since  the  laboratory  targets  where  optically  very  rough, 
the  resultant  contrast  ratio  of  both  speckle  patterns  is  probably  the  same  in 
any  two  experimental  and  simulation  cases  we  wish  to  compare.  The  case 
that  comes  closest  in  other  comparable  parameters  are  the  2 rad/sec  simul- 
ation case  of  Figure  18  and  the  2tt  rad/sec  experimental  case  of  Figure  6. 
Table  4 lists  the  parameters  used  in  each  case.  There  are  a number  of  dif- 
ferences between  the  actual  experiment  and  an  exact  scaling.  First,  the 
simulation  assumes  a uniform  illumination  across  a 10  cm  diameter  spot  at 
the  target.  Inferring  a transmitter  diameter  from  this  gives  = 25  cm, 
which  is  considerably  smaller  than  the  1.  2 m x 1.  0m' annulus  usea  for  the 
receiver.  A fair  amount  of  receiver  averaging  is  therefore  occurring  in  the 
computer  simulation.  On  the  other  hand,  if  a 2.  5 cm  diameter  uniform  spot, 
corresponding  to  a diffraction  limited  1.  0m  transmitter,  were  used  to  illum- 
inate the  target,  the  rotation  rate  would  have  to  increase  by  about  a factor 


of  4,  to  8 rad/ sec,  to  produce  about  the  same  power  within  the  COAT  system 
dither  band.  The  experimental  illumination  function  is  roughly  Gaussian. 


Actually  closer  to  (r)t  but  the  distinction  is  not  important. 


TABLE  4.  EXPERIMENTAL  AND  COMPUTER  SIMULATION  PARAMETERS 
FOR  MOST  CLOSELY  CORRESPONDING  CASE 


Parameter 

Simulation 

Experiment 

wavelength 

10.6  pm 

0.488  pm 

Z,  range 

2 km 

lm 

D^.,  Transmitter 
aperture  diameter 

= 25  dm  (inferred  from 
uniform  illumination 
of  10  cm  spot  at  tar-  , 
get) 

0.  15  mm 

Dd,  Receiver  aperture 
R. 

1 .2  m x 1 .0  m annulus 

0.  15  mm,  filled 
aperture 

= Target  Rotation  Rate 

2 rad/sec 

Ztt  rad/sec 

Target  Surface 

Completely  diffuse 
(randomizes  reflected 
phase  front) 

. ' 

Completely  diffuse 

The  uniformly  illuminated  target  gives  rise  to  a linear  dependence  of  speckle 
spectral  density  on  spatial  frequency  whereas  the  Gaussian-illuminated  tar- 
get produces  a Gaussian  dependence.  The  receiver  aperture  integration  or 
weighting  function  for  the  circular  experimental  aperture  and  the  annular 
simulation  aperture  also  differ. 

The  point  to  be  made  here  is  that  since  the  various  parameters  do  not 
correspond  exactly,  it  is  likely  that  the  target  rotation  rate  required  to 
produce  a given  effect  experimentally  will  differ  from  that  given  by  the  scal- 
ing rules.  In  fact,  this  was  the  case.  Using  the  scaling  rules  in  Eq.  (5) 
directly  gives  3 rad/sec  for  the  experimental  rotation  rate  that  should  give 
speckle  effects  corresponding  to  the  2 rad/sec  simulation  case.  This  rota- 
tion is  roughly  a factor  of  2 lower  than  the  2tt  rad / sec  rate  used.  Neverthe- 
less, when  one  compares  the  two  power  spectra  (see  Figure  6 and  18),  which 
is  really  the  crucial  consideration,  we  see  they  are  almost  identical,  and 
consequently  the  resulting  convergence  levels  are  also  almost  identical.  The 
10.  6 pm  scenario  computer  simulation  shows  an  average  convergence  level 
of  0.  32  with  oscillations  between  0.  22  and  0.42.  The  experimental  observa- 
tion in  Figure  6 has  a mean  of  0.  31  with  oscillations  roughly  between  0.  18 
and  0.  55.  An  additional  observation  that  lends  credence  to  the  comparison 
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is  that  in  both  cases,  increasing  the  rotation  rate  did  not  cause  any  more 
degradation  in  convergence  level.  This  indicates  that  both  power  spectra 
correspond  to  rotation  rates  that  produce  maximum  degradation  in  COAT 
convergence  levels. 

We  have  also  performed  a series  of  experiments  to  observe  the  effect 
of  receiver  size.  In  these  experiments,  the  ratio  of  receiver  to  transmitter 
diameter  was  varied  and  the  normalized  convergence  levels  were  recorded 
with  the  target  rotating  at  2tt  rad/sec.  Table  5 shows  the  results  of  these 
experiments.  Beginning  with  equal  transmitter  and  receiver  apertures, 
doubling  of  the  receiver  aperture  raises  the  convergence  level  from  0,  30  to 
0.  50  (a  factor  of  1.  66).  Referring  to  Figure  A-7  in  the  Appendix,  we  see 
that  Cg  decreases  by  a factor  of  0.  74  for  this  change  in  convergence  level. 
Since  the  speckle  noise  power  is  proportional  to  the  product  of  the  square  of 
Cg  and  the  convergence  level  (see  Eq.  A-27  in  Appendix),  we  can  infer  a 
reduction  of  the  speckle  signal  power  by  a factor  of  0.  66  when  the  receiver 
diameter  is  doubled.  Thus  aperture  integration  does  indeed  help  to  overcome 


speckle  effects. 

TABLE  5 EFFECTS  OF  RECEIVER  APERTURE  SIZE 


SECTION  IV 


CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  WORK 

The  purpose  of  these  investigations  was  to  develop  an  understanding 
of  the  interactive  effects  of  a multidither  COAT  servo  with  the  spurious 
modulations  induced  by  speckle.  We  have  demonstrated  experimentally  and 
verified  theoretically  that  speckle  modulation  generated  by  the  motion  of 
ordinary  target  surfaces  can  under  some  conditions  degrade  the  convergence 
behavior  or  corrective  ability  of  a COAT  system.  Factors  such  as  the  ratio 
of  receiver  size  to  transmitter  size,  target  velocity  or  rotation  rate,  and 
COAT  system  feedback  loop  gain  and  clipping  level  all  influence  the 
susceptibility  of  the  COAT  system  to  speckle. 

Some  of  the  sensitivity  to  scenario  and  COAT  system  parameters  is 
shown  in  the  following  examples:  (1)  Under  speckle  modulation  conditions 

where  the  COAT  system  is  severely  degraded  to  a level  of  0.  3 of  the  maximum 
convergence  level,  doubling  the  receiver  diameter  increases  the  convergence 
level  to  0.  5;  (2)  Under  the  condition  of  high  feedback  loop  gain,  increasing 
the  target  rotation  rate  increases  the  degradation  until  a maximum  is 
reached.  Increasing  the  rotation  rate  beyond  this  point  leads  to  less  degrada- 
tion. There  thus  exists  a rather  narrow  "window"  of  target  motions  that 
produce  significant  problems;  (3)  With  maximum  speckle  degradation,  reduc- 
ing the  COAT  servo  loop  gain  by  10  to  15  dB  usually  removes  most  of  the 
convergence  degradation  but  with  an  increase  in  convergence  time  and 
decrease  in  correction  bandwidth;  Judicious  application  of  clipping  in  the 
servo  loop  in  the  presence  of  speckle  modulation  diminishes  the  convergence 
degradation  with  little  change  in  convergence  time;  and  (4)  Square  wave 
demodulation  and  larger  dither  amplitude  can  also  improve  performance. 

Thus  the  degree  to  which  the  COAT  system  is  affected  by  speckle  is 
very  dependent  on  the  COAT  system  parameters  as  well  as  the  target  param- 
eters. The  crucial  speckle  parameter  Cg,  represents  the  combined  speckle 
power  within  the  servo  bandwidth  around  each  dither  frequency.  The  worst 
case  computer  simulation  results  obtained  correspond  to  Cg  = 0.24  and  an 
average  convergence  level  of  0.  32.  A factor  of  2 reduction  in  speckle 
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contrast  ratio  by  itself  would  reduce  Cg  to  0.  12  and  lead  to  a significant 
increase  in  average  convergence  level  (to  0.6).  In  our  target  signature 
characterizations,  we  found  a factor  of  2 variation  between  the  speckle 
contrast  ratio  of  a metallized  sphere  and  a painted  sphere,  and  a factor  of 
4 variation  in  the  contrast  ratio  of  a metallized  sphere  compared  to  a side- 
illuminated  metallized  cone.  Laser  bandwidth  and  the  ratio  of  target 
roughness  to  optical  wavelength  also  offset  the  contrast  ratio  strongly. 

Another  significant  point  of  interest  is  that  there  was  effectively  no  speckle 
modulation  effect  for  the  nose-on -target  case,  where  the  axis  of  rotation 
coincides  with  the  laser  propagation  direction.  This  is  the  most  probable 
encounter  aspect  in  many  scenarios.  Because  of  these  effects,  it  is  not 
unreasonable  to  expect  less  than  the  theoretical  worst-case  speckle  degrada- 
tion in  a real  target  situation. 

Another  system  consideration  is  that  the  speckle  degradation  is  deter- 
mined by  the  ratio  of  dither  power  to  speckle  power.  By  increasing  the 
dither  power,  we  can  reduce  the  effect  of  speckle.  The  process  is  limited, 
however,  because  convergence  is  also  degraded  by  the  very  large  dither 
amplitudes.  When  other  system  factors  such  as  clipping  level  and  loop 
gain  are  optimally  adjusted  for  a specific  situation,  the  convergence  degrada- 
tion due  to  speckle  modulation  can  also  be  reduced  independently  of  target 

effects.  Another  possible  system  technique  for  reducing  the  effects  of 

7 

speckle  is  optimal  predictive  filtering  (Kalman  filtering). 

As  a final  comment  on  how  system  parameters  can  affect  speckle 
properties,  we  consider  what  happens  when  the  laser  illuminating  the  target 
is  not  monochromatic,  as  is  the  case  for  high-power  HF/DF  (3.8  pm)  and 
CO  (5  (xM)  lasers.  To  examine  the  effect  of  several  laser  lines  illuminating 
the  target  simultaneously,  it  is  necessary  to  know  the  ratio  of  the  rms  sur- 
face roughness,  IT,  to  the  optical  wavelength,  \.  As  shown  below,  target 
surfaces  with  rms  roughnesses  on  the  order  of  0.  1 wave  or  less  can  be 
considered  totally  specular  and  the  speckle  problem  can  be  ignored.  For 
rougher  surfaces,  the  intensity  of  the  speckles  and  their  spatial  frequency 
spectrum  will  have  to  be  calculated. 

Under  certain  circumstances,  the  presence  of  multiple  laser  lines 
can  be  very  beneficial  in  reducing  the  effects  of  speckle.  For  rough  surfaces, 
each  laser  line  generates  a partially  uncorrelated  speckle  pattern  and  when 


the  multiline  speckles  are  combined  the  effect  is  to  wash  out  the  (single- 
line)  speckle  pattern  in  the  same  way  that  completely  incoherent  radiation 
washes  out  speckle.  The  following  first-order  analysis  illustrates  qualita- 
tively what  might  be  expected  for  an  HF/DF  or  CO  laser. 

Consider  a reflecting  surface  with  rms  roughness  F illuminated  by 
monochromatic  light  of  wavelength  X.  It  can  be  shown  that  the  correlation 
coefficient  for  the  separate  speckle  patterns  from  two  monochromatic 
sources  with  a small  wavelength  difference  AX  is 


I = 


(6) 


By  defining  £ = e-1  as  a decorrelation  criteria,  the  wavelength  decorrela- 
tion dependence  is  found  to  be 


^ 4Trh 


(7) 


We  see  that  the  decorrelation  wavelength  difference  is  inversely  proportional 
to  the  surface  roughness.  Figure  22  shows  a plot  of  minimum  change  in 
wavelength  needed  to  decorrelate  a speckle  pattern  as  a function  of  surface 
roughness  calculated  using  Eq.  (7). 

A typical  HF/DF  laser  will  oscillate  simultaneously  on  many  lines,, 

The  exact  lines  and  their  relative  strengths  depend  on  many  factors, 
including  the  resonator  coupling  fraction,  but  typical  values  for  the  maxi- 
mum and  minimum  frequency  spacings  between  lines  are  on  the  order  of 
0.  34  pm  to  0.  17  pm,  respectively.  This  gives  AX/X  values  of  0.  09  to 
0.  0033,  respectively.  Referring  to  Fig.  22,  we  see  that  some  of  the  widely 
spaced  adjacent  lasing  lines  become  uncorrelated  for  reflecting  surfaces 
with  rms  roughness  on  the  order  of  1.  0 wavelength  and  that  as  the  surface 
gets  rougher,  more  of  the.  lines  become  uncorrelated;  finally  at  roughness 
on  the  order  of  10  wavelengths,  all  the  adjacent  lines  speckle  patterns  become 
uncorrelated. 


■ 


For  typical  CO  lasers,  , values  for  AX/X  range  from  0.  002  to  0.  085. 
From  Figure  22,  an  rms  surface  roughness  of  h/X  a 50  is  required  to 
decorrelate  the  closely  spaced  lines,  but  a roughness  of  h/X  =:  1.  0 is 
required  for  the  largest  line-pair  spacing. 

We  thus  can  see  that  the  nature  of  the  target  surface  as  well  as  the 
illuminating  laser  bandwidth  play  important  roles  in  determining  the  severity 
of  speckle  modulation  effects. 

In  the  presence  of  blooming  the  convergence  degradation  due  to 
speckle  modulation  was  approximately  the  same  as  in  the  absence  of  blooming 
(once  the  effect  of  a larger  receiver  is  accounted  for  in  our  experiments). 

When  the  effects  of  blooming,  speckle,  and  artificial  turbulence  were  simul- 
taneously present,  the  COAT  signal -to -noise  ratio  was  marginal  at  all  but 
the  highest  input  power  levels.  At  those,  levels,  we  found  that  the  COAT 
system  produced  an  improvement  of  65%  in  peak  target  irradiance  over  the 
uncorrected  beam.  However,  the  degradation  due  to  speckle  in  the  presence 
of  blooming  and  turbulence  was  42%,,  compared  to  approximately  70%  observed 
in  the  earlier  experiments.  This  is  accounted  for  by  the  receiver  aperture 
being  twice  as  large  as  before,  and  by  the  fact  that  the  system  was  operating 
with  poor  signal -to-noise  ratio  and  was  not  strongly  converged  in  the  absence 
of  speckle.  In  conclusion,  the  effects  of  speckle  modulation  with  and  without 
the  presence  of  blooming  and  turbulence  are  qualitatively  and  quantitatively 
the  same. 

Additional  work  is  warranted  to  define  those  scenarios  that  can  pro- 
duce significant  speckle  modulations  for  high-power  COAT  system  designs 
and  to  assess  the  probability  of  occurrence  for  such  scenarios.  The  ability 
of  advanced  techniques  to  remove  or  reduce  speckle  effects  should  be 
explored  for  those  scenarios  where  a significant  problem  exists.  Among  the 
items  not  addressed  quantitatively  by  this  contract  are  the  effects  of  modu- 
lation index  (particularly  automatic  control),  number  of  channels,  location  of 
dithe?,  frequencies,  and  of  allowing  the  speckle  itself  to  change  as  the  COAT 
convergence  level  changes.  Inclusion  of  the  variation  of  servo  bandwidth 
with  convergence  level  would  increase  the  usefulness  of  the  analysis. 
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APPENDIX  A 


THEORETICAL  ANALYSIS  OF  COAT  SERVO  RESPONSE 
TO  SPECKLE  MODULATION 

A.  MODEL  OF  SPECKLE  EFFECTS  ON  A COAT  SERVO 

The  effect  of  speckle  on  a multidither  COAT  servo  can  be  simply 
understood  by  considering  the  diagram  in  Figure  A-l.  The  intensity  of  the 
coherent  radiation,  IR,  reflected  off  a target  glint  is  averaged  over  a COAT 
receiver  aperture.  The  electrical  signal,  S,  from  the  receiver  to  the  COAT 
servo  is  directly  proportional  to  IR  and,  in  the  absence  of  speckle,  IR  is 
directly  proportional  to  IT>  the  intensity  at  the  glint: 

S = K I„,  (no  speckle)  , (A-l) 

where  K is  the  proportionality  constant.  The  signal,  S,  thus  varies 
instantaneously  (neglecting  the  short  propagation  time)  with  any  change  in 
I,p.  In  general  the  overall  effect  of  speckle  is  to  spatially  modulate  '-he 
intensity  at  a given  plane.  The  spatial  modulation  results  from  the  inter- 
ference of  light  reflected  from  surface  irregularities  on  the  target.  Temporal 
modulation  of  the  intensity  at  a point  results  from  target  motion  which  in 
turn  causes  the  speckle  pattern  to  move.  The  photodetector  signal  S is  thus 
amplitude  modulated  at  rates  that  depend  on  the  speckle  lobe  size  and  the 
receiver  geometry  and  with  a modulation  depth  that  depends  on  both  the 
speckle  contrast  ratio  and  the  receiver  geometry.  We  represent  this  modu- 
lation by  a function  M (t)  which  multiplies  the  signal  S: 

s 

S = KMg(t)Irp(t)  (with  speckle)  . (A-2) 

To  be  physically  meaningful,  M must  always  be  non-negati/e  since 

s 

S is  the  output  of  an  optical  detector.  For  this  study  we  have  limited  our 
attention  to  modulation  functions  with  an  average  value  of  unity  over  a 
typical  COAT  system  convergence  time  (typically  1 to  2 msec).  Since 
speckle  is  a wave -interference  effect  and  consequently  an  energy-conserving 


phenomenon,  multiplying  the  intensity  by  M must  not  alter  the  total  energy 
in  the  reflected  radiation  field.  In  other  words  the  average  energy  incident 
on  the  receiver  should  be  the  same  with  or  without  the  modulation  function 
M . The  intensity  at  (he  receiver  is  no  longer  required  to  be  instantaneously 
proportional  to  the  intensity  on  the  target  glint,  but  the  average  intensity 
at  the  receiver  over  a convergence  time  should  still  be  proportional  to  the 
average  intensity  on  the  glint,  and  the  constant  of  proportionality,  K,  in 
Eq.  (A-2)  should  still  be  the  same.  Mathematically,  this  implies  that  the 
mean  value  or  dc  level  of  M must  be  unity. 

Modulation  functions  with  averages  different  from  unity  correspond 
physically  to  cases  where  the  relative  motion  between  the  target  and  receiver 
are  slow.  In  this  case  very  few  speckle  lobes  sweep  across  the  receiver 
aperture  during  a convergence  time.  Consequently  the  speckle  modulation 
frequencies  are  much  lower  than  the  dither  frequencies  and  will  not  interfere 
with  the  convergence  process.  An  automatic  gain  control  (AGC)  can  be 
employed  in  the  servo  loop  to  remove  these  slow  variations. 

Subject  to  the  above  physically  meaningful  constraints,  we  define  th  : 
otherwise  completely  arbitrary  speckle  modulation  function, 


Mg(t) 


dw  a(w)  e 


-iwt 


(A-3) 


To  aid  in  developing  an  analytical  solution  to  the  speckle  problem,  we 
will  examine  the  nature  of  the  signal  S.  We  first  separate  Irp  into  two 
components, 


s *M  + 1 
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| negligible 
j terms 


M 


+ I. 


D 


(A-4) 


where  refers  to  the  mean  intensity  and  I_  refers  to  the  dither  modula- 
M 'u 

tions.  A mathematical  justification  for  Eq.  (A-4)  is  given  in  the  derivation 
of  the  statistical  model  in  Section  B (See  Eqs.  (A  — 11)  and  (A-12).  Similarly, 
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we  can  separate  M without  approximations  into  a mean  value  or  low  fre- 

s 

quency  component  M^.  and  a high  frequency  component  M^. 


M = M.  , + M 
s Mac 


(A-5 


Combining  Eqs.  (A-2),  (A-4),  and  (A-5),  we  have 


S = + M.J^.  + M Ix , + M I_.  + (higher  order  terms)] 

[MM  M D ac  M ac  D ' 6 J 


(A -6 


The  various  filters  and  synchronous  detectors  in  the  COAT  servo  will  pass 
signals  only  near  the  dither  frequencies.  Thus  we  need  to  consider  only  the 
terms  KM^I^,  containing  the  dither  error  information,  and  con- 

taining the  speckle  induced  interference  at  the  dither  frequencies.  The  rest 
of  the  terms  contribute  little  or  no  signal  in  the  dither  band.  The  last  term, 
^ac^D’  *S  imPorfarit  only  at  very  low  covergence  levels  (<  5%)  which,  as  we 
shall  show,  are  nelow  any  steady-state  levels  for  the  worst  cases  considered 
(see  Section  III  of  this  report  also).  Our  goal  is  to  calculate  the  average 
power  in  both  of  these  signals  and  then  relate  their  power  ratio  to  the  COAT 
system  convergence  level.  To  do  this,  we  must  first  develop  a statistical 
model  of  a multidither  COAT  system. 


STATISTICAL  MODEL  OF  AN  N-CHANNEL  MULTIDITHER  COAT 
SYSTEM 


In  order  to  derive  a good  statistical  model,  we  begin  by  deriving  the 
exact  expression  for  the  radiant  intensity  at  a point  (i.  e.  , glint)  on  a distant 
target  illuminated  by  an  N-channel  segmented  COAT  array.  This  system 
representation  is  perfectly  general  for  the  problem  we  are  considering. 

Consider  the  electric  field  of  N separate  waves  having  the  same 
amplitude  and  wavelength  incident  on  a given  point  (i.  e.  , target  ^lint). 


Assume  the  polarization  of  each  wave  is  adjusted  so  that  they  are  all  plane 
polarized  parallel  to  each  other.  Then  the  total  field,  E^.,  is 


ET(t) 


I 


n = 1 


E 

o 


ir  (t) 

n' 


(A-7 


v/here  E is  the  magnitude  of  each  wave  (assumed  equal  for  each  wave),  and 
Fn(t)  is  the  time-varying  phase  of  the  n wave  at  the  glint. 

The  radiant  intensity,  normalized  to  unity  when  all  = 0 (1.  e.  , 

Eq  = 1/N),  may  be  written  as 


N 


N 


E. 


N 


2 2 exp  i[rn(t)  ■ rm(t)]  • (A-8 


n = 1 m=l 


In  a COAT  multidither  array,  the  N phase  angles  are  physically 


th 


dithered  at  different  frequencies.  For  such  an  array,  the  phase  in  the  n 
channel  is  r (t)  = PR(t)  +4jsinu>nt.  In  this  representation,  Pn(t)  is  a mean 
phase  angle  and  ^ is  the  amplitude  of  the  small  dither  deviation  from  the 
mean  (typically  4j<  tt/6  radians).  Substituting  r,,(t)  = pn(t)  + 4*  8^n  wn*;  into 
Eq.  (A-8)  we  can  expand  the  result  in  Bessel  functions. 
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Since  we  are  only  interested  in  the  modulations  at  the  dither  frequencies  and 
the  overall  convergence  level,  we  need  to  consider  only  the  first  two  terms 
in  the  Bessel  expansions: 


-I  + _I_ 

N N2 
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(A-10) 


Thus  we  can  identify  the  mean  convergence  level,  and  1^,  the  modula 
tions  at  the  dither  frequencies  as: 


M = 


-1  + 

N 2 

In 


J*(+>  N 
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2 2 c°b  <p„  - pm> 


n = l m = l 
n^m 


(A- 11) 


and 


4J0(^)  Jj(ip) 
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N N 


J 2 Sin  ((3n  “ Sinwnt  . (A-12) 


n = l m=  1 
n^m 


1^31 
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We  thus  have 
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■ V$ . 
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I + I + Higher  order  products  of  Bessel  functions  , (A -4) 
M D b 


which  is  identical  to  Eq.  (A-4).  For  small  values  of  i|i,  the  higher  order 
terms  are  negligible. 

The  convergence  level  and  the  dynamics  of  a COAT  system  are 
functions  of  N degrees  of  freedom,  i.  e.  , the  N phase  angles  Pj,  P2>  ^N" 

Given  the  value  of  each  p and  its  time  derivative  at  any  instant  of  time,  we 
can  calculate  the  values  of  any  dynamic  variables  of  interest.  The  average 
or  mean  value,  however,  can  be  calculated  quickly  by  using  the  following 
statistical  model.  Let  a quasi-static  or  steady  state  convergence  level  be 
characterized  by  a random  distribution  of  the  phase  angles  P^  through  p^ 
over  some  interval  a,  as  shown  in  Figure  A-2.  As  a approaches  zero, 
the  COAT  array  will  converge  on  the  target  glint  (i.  e.  , Lj,-*  1 and  thus 
X — *1).  Conversely,  as  a approaches  2 it,  the  array  approaches  a minimum 
convergence  level.  The  range  of  all  possible  relative  phase  differences, 


B - P , between  any  two  channels,  at  a specified  convergence  level,  is 

rn  . _ - 1 


thus  0 to  a.  The  probability  density  that  a given  Pn  - P is  just  P(P)  a 
over  this  range  and  zero  everywhere  else.  The  probability  density  P(Ap), 
for  an  arbitrary  Pm  - pn  = Ap,  is  the  autocorrelation  function  of  P((3), 

given  by: 


P(AP) 


1 Iap 

a 2 

a 


(A-13) 


This  expression  holds  in  the  region  -o'<Ap<a.  Outside  that  region 
P(AP)  = 0.  This  model  allows  us  to  calculate  expectation  values  of  various 
dynamic  variables  of  particular  interest  for  speckle  effects  such  as  the  ^ 
convergence  level,  <IM>,  and  the  rms  of  the  dither  modulations,  < ID  > 
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Using  the  definition  of  in  Eq.  (A-10),  we  write  the  expectation 


value  of  I,,  as: 
M 


T2, . . N N 
JJ+)  v- 


1 " _ \ Y / V—  V~ 

<IM>=N-+-VI  1 <C0S  <Pn  - Pm»  • <A-14> 


n-1  m = 1 
n^m 


Using  Eq.  (A-12),  it  is  easy  to  show  that 


/CO 

P(Ap)  cos  (A p)  d(Ap)  = ~ (1  - cos  a)  . (A-15) 


Thus  we  have 
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. (A-16) 


Figure  A-3  shows  a plot  of  <IM>»  versus  a for  an  18-element  array,  and  a 
peak  dither  amplitude  of  ^ = 0.  349  radians  (20°). 

A similar  calculation  using  Eq.  (A-12)  shows  that  = and 
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(A- 17) 


m = l 
m^n 


In  Eq.  (A- 17)  we  ignore  the  additional  oscillatory  terms  at  frequencies 


u>  + u>  and  to  - ui  because  of  the  assumed  quasi-static  steady  state  where 
n m n m 


the  convergence  level  remains  fairly  constant  over  a long  period  of  time. 
These  terms  will  therefore  have  a time  average  close  to  zero.  Let 


= 2 sin  <p„  - 


(A- 18) 
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Figure  A-3.  Statistical  expectation  value  of 
the  convergence  level  versus  a 
for  an  18-channel  muitidither 
COAT  system  (dither  amplitude 
= 20  ),  calculated  from  eq . 16. 
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Since  b is  a function  of  the  random  variables  (3,,  P?,  . . . > PN>  it  is  also  a 
n 1 ^ 

random  variable  whose  expectation  value  is  expressed  as 
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./  dp>  i 


dpNP(P1)P(P2)  ..  P(PN)  bn(Pj.  P2 


(A-19) 


We  have  assumed  that  the  P's  are  randomly  distributed  between  the  phase 
angles  a and  a + a as  indicated  in  Figure  A-2.  Since  aQ  is  arbitrary,  it 
will  be  convenient  to  define  it  as  a a^  = -a/  2;  this  will  make  the  integrals  in 
Eq.  (A--9)  easier  to  calculate.  Using  the  above  information  plus  the  defini- 
tion of  bn  in  Eq.  (A- 18)  we  can  rewrite  Eq.  (A-19). 
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Carrying  out  the  integrals  gives 
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Substitution  of  Eqs.  (A-18)  and  (A-21)  into  Eq.  (A-17)  gives 


The  rms  dither  amplitude  is  just  the  square  root  of  Eq.  (A-20).  A plot  of 
the  rms  dither  amplitude  versus  a for  an  18 -element  COAT  array  with 
ifj  - o.  349  rad  (20°)  is  shown  in  Figure  A-4.  Since  a is  just  a statistical 
parameter,  it  is  more  meaningful  to  plot  the  dither  amplitude  as  a function 
of  the  convergence  level,  <IM>  as  shown  in  Figure  A-5. 

C.  PREDICTIONS  OF  EQUILIBRIUM  CONVERGENCE  LEVELS 

At  this  point,  we  make  an  assumption  about  the  convergence  of  a 
COAT  system.  Later  we  will  verify  the  assumption  with  computer  simulation 
data.  The  assumption  is  that  the  COAT  array  will  converge  to  an  equilibrium 
level  where  the  average  power  in  the  synchronously  detected  dither  signals, 
PD,  is  equal  to  the  average  power  in  the  synchronously  detected  speckle 
induced  signals,  Pg.  During  the  original  development  of  this  analysis,  we 
chose  to  calculate  the  average  power  in  these  two  signals  before  synchronous 
detection.  It  is  just  as  easy  to  do  the  analysis  this  way  as  long  as  we  recog- 
nize that  the  detection  process  alters  the  relative  power  in  these  signals. 

This  alteration  happens  because  the  phases  of  the  synchronous  detectors  are 
synchronized  to  the  dither  phases  while  the  speckle  induced  modulations  are 
random.  The  detection  process  thus  introduces  a time  averaging  into  the 
speckle  signal  v/hich  does  not  occur  in  the  dither  signal.  The  net  result  is  a 
factor  of  2 increase  in  the  power  ratio. 
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Figure  A-4.  Statistical  expectation  value  of  the  rms 
dither  amplitude  versus  a for  an  18- 
channel  multidither  COAT  system 
(dither  amplitude  = 20°). 


RMS  DITHER  AMPLITUDE  <1 


If  wo  define  the  power  ratio,  p,  before  synchronous  detection 


<PD> 

<PS> 


(A-23) 


then  we  can  define,  p1,  the  power  ratio  after  synchronous  detection  as: 
p1  = Zp.  Thus  the  equilibrium  convergence  level  will  correspond  to  that 
level  where  p = 0.  5.  We  now  proceed  to  calculate  this  predicted  conver- 
gence level. 

Let  SD  be  that  portion  of  the  signal  S containing  the  dither  error 

signal: 


SD  = K MM:D 


(A-24) 


The  expectation  value  of  the  average  power  at  the  dither  frequencies  is 


< PD> 
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where  we  have  used  < > =*  <a(0)^>  = 1.0.  Let  Sg  be  that  portion  of 

the  signal  S in  the  term  KMacI^  which  causes  spurious  modulations  near 
the  dither  frequencies.  In  order  to  calculate  this  signal,  we  define  a 
"speckle  coefficient,  " Cg,  as 


1/2 


Sj 


w.+Au) 
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aHa  H.  du> 
2T 


±u.  to, -Aw 
l l 
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(A-26) 


I his  coefficient  is  the  square  root  of  the  sum  of  2N  integrals,  where  N is 
the  number  of  dither  frequencies  in  the  COAT  servo.  The  summation  is 
carried  out  over  both  positive  and  negative  values  of  the  dither  frequencies, 
wr  The  quantity  Aw  is  the  open  loop,  unity-gain  bandwidth  of  each  COAT 
servo  channel.  These  quantities  are  shown  graphically  in  Figure  A-6  for  the 
positive  half  of  the  power  spectrum  of  an  arbitrary  modulation  function.  The 
shaded  area  corresponds  to  the  integrals  in  Eq.  (A-26). 

Physically,  Cg  corresponds  to  that  portion  of  the  ac  modulation  power 
in  Mg  which  is  at  or  near  enougn  to  each  dither  frequency  so  that  it  can 
create  false  error  signals  which  compete  with  the  correct  dither  signals  for 
control  of  the  convergence  process.  Thus  the  average  power  in  the  signal  Sg 
is  given  by 


= 2‘T  L & 


dt  = K 


-T 


(A-27) 


The  quantity  Cg  has  the  interesting  property  that  it  is  never  larger  than 


0.  707.  This  can  be  understood  by  considering  the  total  average  power,  PTg, 

o'-  the  speckle  modulation  function  M : 

s 


TS 


lim 
T—  oo  2T 


kf 


Mg(t)  Mg(t)  ' dt 


Using  the  definition  of  M (t)  in  Eq.  (A -3)  gives: 


lim 


TS  “ T — * 21 


* /. 


a(u>)  a ’ (u>)  du> 
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(A-29) 
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A limiting  value  for  P_,  occurs  when  M consists  of  a unit  dc  level  and  a 

unit  amplitude  single  frequency  modulation.  In  this  case  PTS  =1.5.  A 

more  realistic  modulation  function  would  include  many  frequencies  with 

random  phases.  The  requirements  that  M must  always  be  non  negative  and 

nave  an  average  value  of  1.  0 however,  prevents  P-j-g  from  exceeding  this 

limiting  value.  Furthermore,  if  the  dc  power  is  always  1.0,  the  ac  power 

will  never  exceed  0.5,  Clearly,  the  sum  of  the  integrals  in  Eq.  (A-26)  will 

always  be  less  than  or  equal  to  the  integral  in  Eq.  (A-29).  In  physical  terms, 

'.his  is  the  obvious  statement  that  the  portion  of  the  ac  speckle  modulation 

power  at  or  near  the  dither  frequencies  is  always  less  than  or  equal  to  the 

total  ac  power.  Since  the  total  ac  power  is  less  than  or  equal  to  0.  5,  then 

C must  also  be  less  than  or  equal  to  0.  5 (and  Cc  5 0.  707). 
s ^ 

We  have  shown  how  to  calculate  Pg  and  PD  as  a function  of  Cg  and 

convergence  level,  < IM>.  It  is  a simple  step,  therefore,  to  calculate  <IM> 

as  a function  of  C for  constant  p.  When  p = 0.  5 the  result  of  this  calcula- 

s 

tion  is  the  analytically  predicted  convergence  level.  In  Figure  A-7  we  show 
the  predicted  level  for  an  18-element  COAT  array  with  dither  phase  modu- 
lation amplitudes  of  = 0,  349  radians  (20  degrees). 


D.  VERIFICATION  OF  ANALYSIS  USING  COMPUTER  SIMULATIONS 


Computer  simulations  have  verified  that  the  analytical  approach 
described  above  gives  very  accurate  predictions  of  average  convergence  levels, 
provided  the  servo  bandwidth  variation  with  convergence  level  can  be  ignored. 

In  a typical  COAT  servo,  the  gain  in  each  channel  is  a function  of  all  the 
phases  in  the  N channels;  i.  e.  , the  gain  is  a function  of  the  overall  system 
convergence  level  and  the  phase  error  in  that  channel.  But  since  the  servo 
bandwidth  is  a function  of  the  gain,  it  also  varies  with  the  overall  convergence 
level  and  the  phase  error  in  that  channel.  Typically,  the  bandwidth  in  all 
channels  is  small  at  low  convergence  levels  and  increases  to  a maximum 
value  at  full  convergence  for  a servo  without  an  AGC.  When  a divider 
AGC  is  included  in  the  loop,  variation  in  gain  and  thus  bandwidth  from 
channel  to  channel  is  very  large  at  low  convergence  levels  and  then  approaches 
a constant,  generally  lower  level  in  all  channels  at  full  convergence.  This  . 
dependence  of  servo  bandwidth  on  convergence  level  makes  a straightforward 
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analytical  prediction  of  convergence  level  somewhat  more  difficult  because 

the  servo  bandwidth  must  be  known  to  calculate  C . The  equilibrium 

s 

convergence  level  is  then  found  as  shown  from  Figure  A-7.  The  servo 
bandwidth  associated  with  this  predicted  level,  however,  may  not  be  the 
same  as  the  bandwidth  originally  chosen  to  calculate  Cs>  Thus,  several 
iterative  calculations  are  necessary  to  insure  that  all  these  parameters  are 
compatible  with  each  other. 

In  order  to  avoid  these  iterative  calculations  and  still  be  able  to 

compare  the  theory  with  computer  simulations,  we  can  alternately  define  an 

artificial  speckle  modulation  function  M as  follows:  . 
r s 


where: 


N = number  of  channels 


M = number  of  erroneous  signals  near  each  dither  frequency 


w = the  jth  frequency  in  a band  of  M equally  incremented 
frequencies  centered  at  the  dither  frequency 

= a randomly  chosen  phase  angle  between  0 and  2tt 


a randomly  chosen  amplitude  between  0 and  an  adjustable 
maximum,  a 

max 


The  to.,  are  chosen  to  lie  in  a sufficiently  narrow  band  around  each 
jk 

dither  frequency  so  that  the  servo  bandwidth  is  always  larger  than  this  band 

at  all  convergence  levels.  A speckle  coefficient  Cg  that  is  effectively 

independent  of  servo  bandwidth  can  then  be  calculated  and  the  resulting 

analytically  predicted  convergence  levels  can  be  compared  to  a computer 

simulation  run  which  uses  an  identical  modulation  function.  In  this  study  we 

chose  21  equally  spaced  frequencies  over  a bandwidth  of  ±83  Hz  around  each 

of  18  dither  frequencies  (i,  e.  , N = 18,  M = 21).  By  choosing  different  values 

for  a , we  were  able  to  gather  data  for  different  values  of  C (C  « 8 a ■ 

max  6 s s max 


Simulations  were  run  at  three  values  of  broadband  signal -to -noise  ratio 
(receiver  shot  noise):  10^,  40,  and  20.  When  (S/N)g.nQt  - 10  and 

C ^ 0 (no  speckle  noise),  the  simulation  produces  a time  history  of  the 
convergence  level,  or  intensity  on  the  target  glint  as  shown  in  Figure  A-8(a). 
When  the  speckle  noise  modulation  shown  in  Figure  A-8(c)  is  present 
(C  = 0.  127),  the  time  history  in  Figure  A-8(b)  results. 

The  behavior  shown  in  Figure  A-8(b)  for  Cg  ^ 0 is  typical.  In 
general,  the  speckle  modulations  have  two  effects:  (1)  the  mean  convergence 

level  is  lowered;  and  (2)  a random  oscillation  of  the  intensity  occurs  between 
a maximum  convergence  level  and  a minimum  level.  This  oscillation  is 
much  slower  than  the  dither  oscillations  (on  the  order  of  0.  5 to  1.  0 kHz 
versus  10.  0 to  40.  0 kHz)  and  its  magnitude  varies  with  convergence  level. 
Figure  A-9  shows  a plot  of  the  maxima  and  minima  recorded  during  these 
simulations  as  a function  of  Cg. 

Also  shown  in  Figure  A-9  is  the  analytically  predicted  mean  conver- 
gence level  corresponding  to  P = 0.  5 and  two  other  analytical  curves  corre- 
sponding to  p = 0.  25  and  p = 1.  0.  Notice  that  the  maximum  and  minimum 
convergence  levels  correspond  very  well  to  these  two  values  of  p.  In  terms 
of  the  power  ratio  after  demodulation,  p',  we  see  that  p'  = 1 corresponds  to 
the  mean  convergence  level,  p'  =0.5  corresponds  to  the  maximum  level, 
and  p'  = 2 corresponds  to  the  minimum  level.  Evidently  when  the  relative 
power  in  the  dither  signal  after  demodulation  is  a factor  of  2 greater  than  the 
speckle  signal,  it  will  totally  dominate  and  the  COAT  system  will  converge, 
conversely,  when  the  power  in  the  speckle  signal  is  a factor  of  two  greater 
than  the  dither  signal  it  will  totally  dominate  and  the  COAT  lyst'em  will  fail 
to  converge  any  further. 


BROAD  BAND  SIGNAL 
TO  NOISE  RATIO  - 101 


NO  SPECKLE 


BROAD  BAND  SIGNAL 
TO  NOISE  RATIO  = 101 


WITH  SPECKLE 
Co  = 0.127 


APPENDIX  B 


SYSTEM  SAFETY  ANALYSIS 

Paragraph  4.4  of  the  contract  statement  of  work  requires  that  a 
system  safety  analysis  be  performed  in  accordance  with  paragraph  5.  8. 2.1 
of  MIL-STD-882.  Such  an  analysis  was  performed  during  the  course  of  the 
contract  work  to  identify,  and  correct  where  necessary,  possible  hazards 
associated  with  the  operation  of  the  system  known  as  the  "multidither  COAT 
experimental  model"  (hereafter  called  "the  COAT  system")  used  in  the 
experimental  tasks  of  this  contract.  The  analysis  performed  follows  a 
similar  analysis*  performed  for  the  same  COAT  system  when  it  was  located 
in  a different  area. 

There  are  two  potential  hazard  sources  for  users  of  this  system 
within  the  laboratory  environment  where  this  contract  work  was  performed. 

(1)  electrical  shock  and  (2)  eye  damage  by  laser  radiation.  Appropriate 
safety  precautions  have  been  taken  as  part  of  the  Hughes  Industrial  Safety 
Program  to  insure  sufficient  fire  protection,  clearance  for  egress  from 
the  laboratory,  access  to  power  line  breakers,  etc.,  so  that  the  require- 
ments of  the  California  OSHA  codes  are  satisfied. 

Electrical  shod  'om  the  COAT  system  was  assessed  to  be  in 
Category  I - Negligible  as  defined  in  paragraph  3.14  of  MIL-STD-882. 

Proper  design  of  the  electrical  parts  of  the  system  insures  that  no  personnel 
can  encounter  lethal  or  hazardous  electrical  voltage  and/or  current  levels. 

The  haza.rd  level  for  optical  eye  damage  was  assessed  to  be  in 
Category  II  — Marginal.  The  nature  of  the  program  required  working  with 
two  laser  sources,  a 30  mW  He-Ne  unit  operating  at  0.6328  pm  and  a 1.5  W 
argon  unit  with  an  output  wavelength  of  0.488  pm.  Various  system  optics 
reduce  the  maximum  potential  power  levels  to  10  mW  at  0.6328  pm  and 
100  mW  at  0.488  pm.  The  potential  hazard  is  adequately  controlled  by  pro- 
viding safety  glasses  for  all  personnel  while  working  on  or  around  the  system 
and  by  using  masks,  optical  stops,  etc.  to  confine  the  optical  laser  beams  to 


*W.B.  Bridges,  etal.,  RADC-TR-74-108,  Jan.  1974,  Contract 
F30602-73-C-0248,  (779668). 


areas  where  no  personnel  can  directly  view  them  at  day  point  where  the 
optical  power  density  exceeds  2.5  mW/cm  . 

Sufficient  safety  procedures  have  thus  been  implemented  to  insure 
that  both  the  0.6328  pm  radiation  (a  Class  III  hazard  at  the  laser  output)  and 
the  0.488  pm  radiation  (a  Class  IV  hazard  at  the  laser  output)  are  adequately 
controlled  to  reduce  the  hazard  to  Category  II  for  personnel  working  within 
the  laboratory  where  the  COAT  system  is  operated.  In  addition,  personnel 
are  instructed  in  proper  safety  procedures  while  within  the  laboratory. 
Appropriate  warning  signs  notify  visitors  to  the  laboratory  and  laser  safety 
goggles  are  available  for  visitors. 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research , exploratory  and  advanced 
development  programs  in  command,  control,  and  communications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability , maintainability  and 
compatibility. 
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